{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Jacobi"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "include(\"../src/Jacobi.jl\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "dfun (generic function with 1 method)"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#fun(x,k=1) = cos(k*x) + 1\n",
    "#ifun(x,k=1) = x + sin(k*x) / k\n",
    "fun(x) = cos(x)\n",
    "ifun(x) = -x.*sin(x) + sin(x) - cos(x)\n",
    "dfun(x) = -sin(x)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "1.682941969615793013305004643260597999245126121596742131345503419983820808782497"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "function calc_integr{T<:Number, QT<:Jacobi.QUADRATURE_TYPE}(Q, a, b, ::Type{QT}=Jacobi.GLJ, ::Type{T}=Float64)\n",
    "    \n",
    "    z = Jacobi.qzeros(QT, Q, a, b, T)\n",
    "    w = Jacobi.qweights(QT, z, a, b)\n",
    "    \n",
    "    y = fun(z)\n",
    "    \n",
    "    return sum(w .* y)\n",
    "end\n",
    "\n",
    "calc_integr(50, 1, 0, Jacobi.GLJ, BigFloat)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "1.682941969615793013305004643260597999245126121596742131345503419983820808782479337897279487086105391708698075815841348586518237841983797762386820655458424818961583911653533213999815528023956817565465132694969605740597312314035924910789787145849340254172972562106764061122754436407736899335523348532478026765491"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "set_bigfloat_precision(1024)\n",
    "\n",
    "Iex{T<:Number}(::Type{T}=Float64) = ifun(one(T)) - ifun(-one(T))\n",
    "Iex(BigFloat)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "Q = 2:40\n",
    "tipo =BigFloat\n",
    "Ie = Iex(tipo)\n",
    "QT = Jacobi.GRJP\n",
    "a = 1\n",
    "b = 0\n",
    "err = [abs(calc_integr(q, a, b, QT, tipo) - Ie) for q in Q];"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "using PyPlot"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAr8AAAIQCAYAAACIUwbVAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3X9U1Ped7/EnAzY70SGpOonbaa1mEskg2qpXkNDV3W0wplQjpgpImzJbt7nm3pOa7mVcJY2axnQ79G5Z29veFbNQ6+yIcvSWxB9Rs72r3UyHRNNWHZv0UooNe1YQ3PIjScMw3D9ADMr4gwzznWFej3M4Md/Pd5j3zOl2X/34+b7fSX19fX2IiIiIiCQAk9EFiIiIiIhEi8KviIiIiCQMhV8RERERSRgKvyIiIiKSMBR+RURERCRhKPyKiIiISMJQ+BURERGRhKHwKyIiIiIJQ+FXRERERBKGwq+IiIiIJAyFXxERERFJGGMu/BYWFjJlyhRSU1NJS0ujsrLS6JJEREREJEYk9fX19RldRCSdO3eO++67j5SUFF577TUWLlxIIBBg+vTpRpcmIiIiIgZLMbqASHM4HIN/Tk5OJjU1FYvFYmBFIiIiIhIrxtyxB4Di4mLMZjM5OTls376dyZMnG12SiIiIiMSAMXfs4bJQKERdXR1/9Vd/xS9+8QumTp16zT0XL17k5ZdfZtq0aZjNZgOqFBEREZHreffdd/nd737HQw89FJkNzb44tmvXrr4JEyb0TZgwoe9zn/vcsPd8/vOf7/ve974X9vWAfvSjH/3oRz/60Y9+Yvxn165dEcmPUT3z29XVxbPPPssvfvEL3njjDdra2ti0aRObNm0a9t6nn36avXv30t7ezv3338/f/u3fUlBQMHhPcXExxcXF133PYDDIhAkThl2bNm0aALt27RpyVni0fOELT9LY+A9A0jCrIaZPX0dt7bZRryMSnnrqKb773e8aXUZc0Xc2Mvrebp2+s5HR93br9J2NjL63W3Pu3Dm++MUvDua2Dyuq4ffixYtUVlby6U9/mvz8fHbs2EFS0nBBEFasWMHrr7/Ot7/9bWbMmIHH46GoqIhQKERRUdGwr7lw4QInTpzg4Ycf5rbbbmPfvn34/X527Ngx7P2Xjzo4HA7mzp0bmQ95HePGTQTmXXc9GnVEwh133BE3tcYKfWcjo+/t1uk7Gxl9b7dO39nI6HsbmUgdUY1q+J02bRqXLl0CoK2tLWwoPXjwIMeOHcPr9Q7u9C5atIimpiZKS0spKCjAZBr+Wb1t27axZs0aTCYTGRkZ1NXVYbPZRucD3aKUlF76d+6H3/l9++1ePB5YsQJ0BFlEREQk8gxrddZ3nefs9u/fj8ViYeXKlUOuO51OVq9ejd/vJzs7+5rX3X333Rw/fvyWa8nNzWXcuHHYbLbBoFxUVBR2h3mkMjPTCQT8wIJhVv1YLOl88Ytw553wxS/C00/D3XcPvau1tRWXy019fYBgMJmUlF4yM9Nxu11YrdaI1isiIiISTV6vF6/XC0BzczPNzc309PRE9D1iss/vmTNncDgc1+zuzpo1C4CzZ88OG35H6ujRo1H56we328WJEwU0NGwFsujvNBcC/NjtZfh8Nfznf8I//RPs2QPf/ObQ17e0tPDAA4U0NDwPuOnfQQ4RCNRz4kQBPl+NArCIiIjEreE2H0+dOsW8eeGPjd6qmOzz29bWxsSJE6+5fvlaW1tbtEuKCKvVis9XQ0nJPtLTlzJjxjLS05dSUrJvMLjedx9861vw//5f/w7wB61fXz4QfBdw5eiECVhAQ8NWXC531D5LpHfFE4G+s5HR93br9J2NjL63W6fvbGT0vRkrJnd+xzKr1UpVVfkN7xvuOcD6+gD9O77DyaK+/rkPVdut0P/h3jp9ZyOj7+3W6TsbGX1vt07f2cjoezNWTO78Tpo0adjd3fb29sH1RBQMJjP8w3IApoF1EREREQknJsPv7NmzOXfuHKFQaMj106dPA5CRkWFEWYa70i1iOCGSknqjWY6IiIhI3InJ8Jufn09XVxe1tbVDrldXV2Oz2cjKyjKoMmNlZqYD/jCrfhoa0nn6aRjYIBcRERGRq0T9zO+hQ4fo7u6ms7MT6O/ccDnk5uXlYTabWbJkCbm5uaxdu5aOjg7sdjter5cjR47g8XjCDsYY667XLWLatDLy8mr47ndh2zbYuROWLze2XhEREZFYk9R3vYa7o2D69Ok0NTX1v3lS0mC/36SkJBobG5k6dSoA3d3dlJWVsWfPHtrb23E4HGzYsIFVq1ZFrJbLrTNOnjwZN5NWbtTnt6UFvvMdWLsWpk+/tdeKiIiIxJpI57Woh99YEo/hd6SG9gjO4nKPYKjHbt+oHsEiIiISkyKd12LyzK9EXiz1CBYRERExisJvgujvERzuQcGsgXURERGRsU3hN0GoR7CIiIiIwm/CuFGP4Hff7SVxT3+LiIhIolD4TRA36hH8+9+n8xd/Aa+/Hs2qRERERKJL4TdBuN0u7PaNgI/+Lg8M/NOH3V6G1+vi4kWYPx8COv4rIiIiY1TUh1yIMaxWKz5fzUCf3+eu6vPb3+bsC1+AI0cgPd3oakVERERGh8JvArFarVRVlYddT0mBz30u/Os1JENERETincKv3JShQzLcXB6SEQjUc+JEgYZkiIiISFzQmV+5KRqSISIiImOBwq/cFA3JEBERkbFA4VduioZkiIiIyFig8Cs35UZDMt57rzea5YiIiIiMiMKv3JQbDck4fz6dr389mhWJiIiI3DqFX7kpNxqS8YMfuHjoIePqExEREbkZcRd+f/jDHzJ37lw+8pGPsGXLliFrhYWFTJkyhdTUVNLS0qisrDSoyrHn8pCMkpJ9pKcvZcaMZaSnL6WkZB8+Xw1r11oVfkVERCTmxV2f34997GM8++yz7Ny5k6SkoQ9gbdq0iV27dpGSksJrr73GwoULefDBB5k+fbpB1Y4tNxqScT0akCEiIiKxIO7C7yOPPALAT37yE/r6hj6A5XA4Bv+cnJxMamoqFoslqvXJtTQgQ0RERGJF3B17uJHi4mLMZjM5OTls376dyZMnG11SwtOADBEREYkVYy78ejweuru78Xq9OJ1Ozp8/b3RJCU8DMkRERCRWxHT49Xg8WCwWLBYLeXl5N/06k8nE8uXLycnJoa6ubhQrlJuhARkiIiISKyIafru6unC5XCxevBir1YrJZLqmI8MH7123bh02mw2z2cycOXOoqakZck9xcTGdnZ10dnZy4MCBa37H1Q+8XS0YDDJhwoSRfyCJiBsNyHj/fQ3IEBERkeiIaPi9ePEilZWV9PT0kJ+fD4QPqCtWrGDnzp1s3ryZw4cPM3/+fIqKivB6vdd9j97eXt577z2CwSA9PT289957hEIhLly4QG1tLd3d3QSDQfbs2YPf7yc3NzeSH1FG4EYDMlpa0nn//WhWJCIiIokqot0epk2bxqVLlwBoa2tjx44dw9538OBBjh07htfrpaCgAIBFixbR1NREaWkpBQUFmEzD5/JvfvObPPvss4P/vnXrVqqrq3nooYfYtm0ba9aswWQykZGRQV1dHTabLZIfUUbA7XZx4kQBDQ1b6T/7a6J/QIafe+4pY/fuGj7yEWNrFBERkcQwaq3Orm5D9kH79+/HYrGwcuXKIdedTierV6/G7/eTnZ097Gs3b97M5s2bh107fvz4iGrNzc1l3Lhx2Gy2wbBcVFREUVHRiH6fDHV5QEZ/n9/nrurzqzZnIiIi0s/r9Q6eAmhubqa5uZmenp6IvochfX7PnDmDw+G4Znd31qxZAJw9ezZs+B0NR48eZe7cuVF7v0T0YQZkiIiISGIYbvPx1KlTzJs3L2LvYUj4bWtr4957773m+sSJEwfXRQD6+uDMGZgyRRPiRERE5MOLuwlvklhefBEeeaSF1NRCOjo0IU5EREQ+HEP6/E6aNGnY3d329vbBdRGAvDzIzCwfCL6aECciIiIfjiHhd/bs2Zw7d45QKDTk+unTpwHIyMgwoiyJQcnJ0NWlCXEiIiISGYaE3/z8fLq6uqitrR1yvbq6GpvNRlZWuKAjiUgT4kRERCRSIn7m99ChQ3R3d9PZ2Qn0d264HHLz8vIwm80sWbKE3Nxc1q5dS0dHB3a7Ha/Xy5EjR/B4PDec3CaJ5cqEuOH+cxEaWBcRERG5sYiH3yeeeIKmpiagf7rb3r172bt3L0lJSTQ2NjJ16lQA9u3bR1lZGc888wzt7e04HA52797NqlWrIl2SxLnMzHQCAT/9Z36v5h+YICciIiJyYxEPv42NjTd13/jx46moqKCioiLSJcgYc70JcXZ7GW53jbEFioiISNxQqzOJeZoQJyIiIpGi8Ctx4VYnxPX1QWMj3HPPKBYlIiIicceQbg8io+2llyAtDZ5+Gt5+uxWns5SZM/NIS1vGzJl5OJ2ltLa2Gl2miIiIRJl2fmVMeugh+MY34JvfbMHtLqSnR9PhRERERDu/MkZ95CPwzDOQl1c+EHw1HU5EREQUfmWM+81vNB1ORERErlD4lTFN0+FERETkgxR+ZUy7Mh1uOJoOJyIikmgUfmVM65/+5g+zqulwIiIiiUbhV8Y0t9uF3b4R8NE/FY6Bf/oGpsO5jCtOREREok6tzmRMG8l0uN5eSNZRYBERkTFJ4VfGvFuZDrd/P2zdCj/6EcycOcqFiYiISNTp2IPIB3z849DdDXPnQnl5/y5wa6smxImIiIwV2vkV+YD58+HUqf7pcOvXw549LbS0FHL+vCbEiYiIjAXa+RW5itkM3/kOHD8Ov/51+UDw1YQ4ERGRsSCuwu/777+P0+lk6tSp3HHHHWRnZ+Pz+QbXCwsLmTJlCqmpqaSlpVFZWWlgtRLvPvMZ+PjHNSFORERkLImr8BsMBrnnnnt49dVX+cMf/sDatWtZtmwZ7733HgCbNm3i7bffpqOjg127dvHkk0/S2NhocNUSz0IhTYgTEREZS+Iq/N5+++184xvf4OMf/zgAjz32GKFQiN/85jcAOBwOUlL6jzEnJyeTmpqKxWIxrF6Jf5oQJyIiMrbEVfi92q9//Wveffdd7Hb74LXi4mLMZjM5OTls376dyZMnG1ihxDtNiBMRERlb4jb8vvPOO3zpS1/iG9/4BrfffvvgdY/HQ3d3N16vF6fTyfnz5w2sUuKdJsSJiIiMLTEdfj0eDxaLBYvFQl5e3uD1np4eVq5cSUZGBhs2bLjmdSaTieXLl5OTk0NdXV00S5Yx5vKEuJKSfaSnL2XGjGWkpy+lpGTfkDZn//qvcPKkwcWKiIjIDUU0/HZ1deFyuVi8eDFWqxWTycSWLVvC3rtu3TpsNhtms5k5c+ZQU1Mz5J7i4mI6Ozvp7OzkwIEDAIRCIb70pS/xkY98hBdeeOG69QSDQSZMmBCZDycJ6/KEuLNnD/Dmm3WcPXuAqqryIf19/+f/hAUL4LnnIBg0sFgRERG5roiG34sXL1JZWUlPTw/5+fkAJCUN/6T8ihUr2LlzJ5s3b+bw4cPMnz+foqIivF7vdd/j8ccf58KFC+zevRuT6Ur5Fy5coLa2lu7uboLBIHv27MHv95Obmxu5DygSRm1t/1CMTZtg4UJoaDC6IhERERlORCe8TZs2jUuXLgHQ1tbGjh07hr3v4MGDHDt2DK/XS0FBAQCLFi2iqamJ0tJSCgoKhgTby5qamnjhhRcwm81DHmQ7fPgw9957L9u2bWPNmjWYTCYyMjKoq6vDZrNF8iOKDOsjH+nf9X34YXjsMfjUp6CiApYta2X9ejf19QGCwWRSUnrJzEzH7XZpMpyIiIgBRm28cV9fuPZQsH//fiwWCytXrhxy3el0snr1avx+P9nZ2de87pOf/CShUOia65cdP358RLXm5uYybtw4bDbbYFguKiqiqKhoRL9PEldODvziF7BuHfz1X7fwta8V8s47Go0sIiJyM7xe7+ApgObmZpqbm+np6Ynoe4xa+L2eM2fO4HA4rtndnTVrFgBnz54dNvyOlqNHjzJ37tyovZ+MbRYLvPACNDaW89OfXh6NfNnQ0chVVeUGVSkiIhJ7htt8PHXqFPPmzYvYexjS7aGtrY2JEydec/3ytba2tmiXJBJxFy5oNLKIiEisielWZyLxrH/0sUYji4iIxBJDwu+kSZOG3d1tb28fXBeJdxqNLCIiEnsMCb+zZ8/m3Llz1zy8dvr0aQAyMjKMKEskom5mNPJ1ngsVERGRUWBI+M3Pz6erq4va2toh16urq7HZbGRlhTsnKRI/bjQa+emnXWRlwYsvGlejiIhIool4t4dDhw7R3d1NZ2cn0N+54XLIzcvLw2w2s2TJEnJzc1m7di0dHR3Y7Xa8Xi9HjhzB4/GEHYwhEk8uj0Z2udzU1z93VZ/fGkwmK1OmwLJl8MQT8J3vgNlsdNUiIiJjW1Lf9RryjsD06dNpamrq/+VJSYP9fpOSkmhsbGTq1KkAdHd3U1ZWxp49e2hvb8fhcLBhwwZWrVoVyXKu63LrjJMnT6rVmRiirw/+9/+Gr38d7rkHvF6YPdvoqkRERGJHpPNaxI89NDY2EgqFCIVC9Pb2Dvnz5eALMH78eCoqKvj3f/933nvvPd54442oBl+RWJCUBGvXwuuvQ0oKZGbCtm3oLLCIiMgoUaszkRgwcyb4/fBf/yt87WuQlwcdHdDa2orTWcrMmXmkpS1j5sw8nM5SWltbjS5ZREQkLhky4U1ErvUnfwIVFfDQQ7BzJ7zzTguf+UwhDQ0ajywiIhIp2vkViTEPP9x/9nfDhvKB4LuAK8Myho5HFhERkVuj8CsSo/rHH2s8soiISCQp/IrEKI1HFhERiTyFX5EYpfHIIiIikafwKxKjbjQeefr09GiWIyIiMiYo/IrEqOuNR77ttjIOH3ZRXg6hUPjfISIiIkMp/IrEqMvjkUtK9pGevpQZM5aRnr6UkpJ9/Pa3NfzN31hxufp7Are0GF2tiIhIfFCfX5EYZrVaqaoqH3bt29+Gv/xLeOwx+NSnYNcu+Oxno1ygiIhInNHOr0gce+gh+OUvISMDcnNh+3ajKxIREYlt2vkViXNTpsDLL8N3vqOdXxERkRtR+BUZA0wmcLmGXmttbcXlclNfHyAYTCYlpZfMzHTcbpfGIouISMJS+BUZg1paWnjggcKB8chu+odlhAgE6jlxogCfr0YBWEREEpLO/IqMQevXlw8E3wVcmRJnAhbQ0LAVl8ttXHEiIiIGirvwW1hYyJQpU0hNTSUtLY3KysqbWhNJJPX1ASArzGrWwLqIiEjiibvwu2nTJt5++206OjrYtWsXTz75JI2NjTdcE0kkwWAyV3Z8r2YaWBcREUk8cXfm1+FwDP45OTmZ1NRULBbLDddEEklKSi/Qx/ABOERvb2+UKxIREYkNcbfzC1BcXIzZbCYnJ4ft27czefLkm1oTSRSZmemAP8yqn9/9Lp0f/SiaFYmIiMSGuAy/Ho+H7u5uvF4vTqeT8+fP39SaSKJwu13Y7RsBHxAauBoCfEyfXsaqVS5KSsDphO5uw8oUERGJupgOvx6PB4vFgsViIS8vb8iayWRi+fLl5OTkUFdXd9NrIonAarXi89VQUrKP9PSlzJixjPT0pZSU7MPvr+Gf/9nKj34Ee/bAwoWgUxAiIpIoIhp+u7q6cLlcLF68GKvVislkYsuWLWHvXbduHTabDbPZzJw5c6ipqRlyT3FxMZ2dnXR2dnLgwIFhf08wGGTChAm3vCYy1lmtVqqqyjl79gBvvlnH2bMHqKoqH+zv+9hj8PrrsH49JOv5NxERSRARDb8XL16ksrKSnp4e8vPzAUhKGv6J8xUrVrBz5042b97M4cOHmT9/PkVFRXi93rC//8KFC9TW1tLd3U0wGGTPnj34/X5yc3OvuyYiw3M4YNUqo6sQERGJnoh2e5g2bRqXLl0CoK2tjR07dgx738GDBzl27Bher5eCggIAFi1aRFNTE6WlpRQUFGAyDZ/Lt23bxpo1azCZTGRkZFBXV4fNZuPChQth10Tk1mk8soiIjEWj1uqsr68v7Nr+/fuxWCysXLlyyHWn08nq1avx+/1kZ2df87q7776b48ePD/s7r7d2I7m5uYwbNw6bzTYYlouKiigqKhrR7xOJdxqPLCIiRvB6vYOnAJqbm2lubqanpyei72FIn98zZ87gcDiu2d2dNWsWAGfPnh02/I6Wo0ePMnfu3Ki9n0isGzoe+bKh45GrqsoNqk5ERMaq4TYfT506xbx58yL2HoZ0e2hra2PixInXXL98ra2tLdolicgHaDyyiIiMVTHd6kxEjKHxyCIiMlYZEn4nTZo07O5ue3v74LqIGOfKeOThhAbWRURE4o8h4Xf27NmcO3eOUCg05Prp06cByMjIMKIsERlwo/HIDkd6NMsRERGJGEPCb35+Pl1dXdTW1g65Xl1djc1mIysr3FlDEYmG641HTkkp48gRF+fOGVefiIjISEW828OhQ4fo7u6ms7MT6O/ccDnk5uXlYTabWbJkCbm5uaxdu5aOjg7sdjter5cjR47g8XjCDsYQkei4PB65v8/vc0P6/JaV1bB7t5UZM4yuUkRE5NYl9V2vIe8ITJ8+naampv5fnpQ02O83KSmJxsZGpk6dCkB3dzdlZWXs2bOH9vZ2HA4HGzZsYFUUx01dbp1x8uRJtToTERERiUGRzmsR3/ltbGy8qfvGjx9PRUUFFRUVkS5BRERERGRYanUmIiIiIgnDkAlvIjJ2vfsunD/fyt/9nZv6+sCQ88Jut0tjkUVExFAKvyISUU891cKOHYX09j4PuOkflhEiEKjnxIkCfL4aBWARETGMjj2ISES1t5cPBN8FXJkSZwIW0NCwFZfLbVxxIiKS8BR+RSSizp4NAOF6dWdRXx+IZjkiIiJDKPyKSEQFg8lc2fG9mmlgXURExBgKvyISUSkpvUC49uGhgXURERFjKPyKSERlZqYD/jCr/oF1ERERYyj8ikhEud0u7PaNgA8IDVwNAT7s9jLcbpdxxYmISMJTqzMRiSir1YrPV4PL5aa+/rmr+vyqzZmIiBhL4VdEIs5qtVJVVW50GSIiItfQsQcRiQnBoNEViIhIIlD4FRHDvfwyfOpTcO6c0ZWIiMhYp/ArIoabNg2SkmD+fNizB1pbW3E6S5k5M4+0tGXMnJmH01lKa2ur0aWKiEic05lfETFcWhr8/Ofw1a9CQUELd9xRyB/+8Dzgpn9gRohAoJ4TJwrw+fTQnIiIjFzc7vz6fD5MJhNbt24dvFZYWMiUKVNITU0lLS2NyspKAysUkVsxYQJ4PJCZWT4QfBdwZVKcCVhAQ8NWXC63cUWKiEjci8vwGwqFeOqpp8jOziYp6coY1U2bNvH222/T0dHBrl27ePLJJ2lsbDSwUhG5FUlJ0NUVALLC3JFFfX0gmiWJiMgYE5fHHv7xH/+RnJwc2tvb6eu7MkbV4XAM/jk5OZnU1FQsFosRJYrICAWDyVzZ8b2aaWBdRERkZOJu57etrY1t27axadOmYdeLi4sxm83k5OSwfft2Jk+eHOUKReTDSEnpBfrCrIYG1kVEREYm7sLvhg0b+Ju/+RtSU1MBhhx7APB4PHR3d+P1enE6nZw/f96IMkVkhDIz0wF/mFX/wLqIiMjIxHT49Xg8WCwWLBYLeXl5nDx5kjfeeIOvfOUrAPT19Q059nCZyWRi+fLl5OTkUFdXF+2yReRDcLtd2O0bAR8QGrgaAnzY7WW43S7jihMRkbgX0fDb1dWFy+Vi8eLFWK1WTCYTW7ZsCXvvunXrsNlsmM1m5syZQ01NzZB7iouL6ezspLOzkwMHDvCzn/2MQCDAXXfdhdVqpaamhm9961uUlJQM+x7BYJAJEyZE8iOKyCizWq34fDWUlOwjPX0pM2YsIz19KSUl+9TmTEREPrSIPvB28eJFKisr+fSnP01+fj47duy45ljCZStWrOD111/n29/+NjNmzMDj8VBUVEQoFKKoqGjY16xZs4aVK1cC/bu+X//615k+fTrr16/nwoULnDhxgocffpjbbruNffv24ff72bFjRyQ/oohEgdVqpaqq/Lr3tLfDM8/A88/DwCkoERGRG4po+J02bRqXLl0C+h9MCxc8Dx48yLFjx/B6vRQUFACwaNEimpqaKC0tpaCgAJPp2k3p8ePHM378+MF/v/3220lNTeWjH/0oFy5cYNu2baxZswaTyURGRgZ1dXXYbLZIfkQRiRG//jX8+Mdw7Bjs3w8faPYiIiIS1qi1OhvuLO5l+/fvx2KxDO7iXuZ0Olm9ejV+v5/s7OwbvkdVVdXgn++++26OHz8+olpzc3MZN24cNpttMCwXFRWF3YEWEeM98AC8/jrk50NmJvzoR7BihdFViYjIh+H1evF6vQA0NzfT3NxMT09PRN/DkD6/Z86cweFwXLO7O2vWLADOnj17U+E3Uo4ePcrcuXOj9n4iEhn33dc/FvkrX4FHH4W//Vv42tda2bDBTX19gGAwmZSUXjIz03G7XTovLCIS44bbfDx16hTz5s2L2HsYEn7b2tq49957r7k+ceLEwXURkZsxYQLs3g3z54PL1cI//EMh7777POCmf1hGiECgnhMnCvTAnIiIxHarMxGRm5GUBP/jf8DixeUDwXcBV6bEmYAFNDRsxeVyG1ekiIjEBEPC76RJk4bd3W1vbx9cFxG5Vb//fQDICrOaRX19IJrliIhIDDIk/M6ePZtz584RCoWGXD99+jQAGRkZRpQlInEuGEzmyo7v1UwD6yIiksgMCb/5+fl0dXVRW1s75Hp1dTU2m42srHA7NyIi4aWk9ALhOs2EBtZFRCSRRfyBt0OHDtHd3U1nZyfQ37nhcsjNy8vDbDazZMkScnNzWbt2LR0dHdjtdrxeL0eOHMHj8YQdjCEicj2ZmekEAn76z/xezU9mZnq0SxIRkRiT1He9hrwjMH36dJqamvp/eVLSYL/fpKQkGhsbmTp1KgDd3d2UlZWxZ88e2tvbcTgcbNiwgVWrVkWynOu63Drj5MmTanUmMga0traSnV1AQ8NW+s/+moAQ4MduL8Pnq6Gx0UpmprF1iojIzYt0Xov4zm9jY+NN3Td+/HgqKiqoqKhZ8ysPAAAgAElEQVSIdAkikqCsVis+Xw0ul5v6+ueu6vNbw7lzVhYtgnXroLwcUgxp9igiIkbSf/WLyJhitVqpqiofdm3yZPj+9/vD7y9/CTU1oLa/IiKJRX1+RSRhJCXBf/tv8MorcOYM/Jf/AqdOGV2ViIhEk8KviCSchQvh5Mn+Xd+cHPB4jK5IRESiReFXRBLSJz4BJ07AqlXwxS/CM88YXZGIiESDzvyKSMIym6G6GubNg+nT+6+1trYOPDAXuOqBORdWHRAWEYl7Cr8iktCSkuDJJ/v/3NLSwgMPFNLQ8Dzgpn9aXIhAoJ4TJwrw+WoUgEVE4pyOPYiIDFi/vnwg+C7gyphkE7CAhoatuFxu44oTEZGIUPgVERlQXx+gfzjGcLIG1kVEJJ4p/IqIDAgGk7my43s108C6iIjEM4VfEZEBKSm9QLiJ76GBdRERiWcKvyIiAzIz0wF/mFU/f/qn6dEsR0RERoHCr4jIALfbhd2+EfABoYGrIcDH7beX8corLr71LegLtzksIiIxT63OREQGWK1WfL6agT6/zw3p8/t3f1fDD35gZeNG+OUv4YUXYPx4oysWEZFbpfArIvIBVquVqqryYde2bIFZs+DLX4bPfAb+z/+BT34yygWKiMiHomMPIiK34AtfAJ8P/vM/oXz4jCwiIjEs7sLvn//5n2M2m7FYLFgsFh588MHBtcLCQqZMmUJqaippaWlUVlYaWKmIjFWzZ8Nrryn8iojEo7gLv0lJSbzwwgt0dnbS2dnJsWPHBtc2bdrE22+/TUdHB7t27eLJJ5+ksbHRwGpFZKyaPBnM5qHXWltbcTpLmTkzj7S0ZcycmYfTWUpra6sxRYqIyDXi8sxvX5hHrR0Ox+Cfk5OTSU1NxWKxRKssEUlgLS0tPPBA4cB4ZDf9wzJCBAL1nDhRgM9Xg9VqNbhKERGJu51fgKeeeoq77rqLBx98kF/+8pdD1oqLizGbzeTk5LB9+3YmT55sUJUikkjWry8fCL4LuDIlzgQsoKFhKy6X27jiRERkUNyFX7fbze9+9zt+//vfk5eXx8MPP0xHR8fgusfjobu7G6/Xi9Pp5Pz58wZWKyKJor4+AGSFWc0aWBcREaPFdPj1eDyDD7bl5eUBMH/+fG6//XZuu+02nnrqKSZPnsyrr7465HUmk4nly5eTk5NDXV2dEaWLSIIJBpO5suN7NdPAuoiIGC2i4berqwuXy8XixYuxWq2YTCa2bNkS9t5169Zhs9kwm83MmTOHmpqaIfcUFxcPPth24MCBW64nGAwyYcKEEX0WEZFbkZLSC4Qb/RYaWBcREaNFNPxevHiRyspKenp6yM/PB/q7MwxnxYoV7Ny5k82bN3P48GHmz59PUVERXq837O//wx/+wNGjR/njH//I+++/z7Zt22hpaSE7O5sLFy5QW1tLd3c3wWCQPXv24Pf7yc3NjeRHFBEZVmZmOuAPs+rH4UiPZjkiIhJGRLs9TJs2jUuXLgHQ1tbGjh07hr3v4MGDHDt2DK/XS0FBAQCLFi2iqamJ0tJSCgoKMJmuzeU9PT1s3LiRN998k3HjxvHpT3+agwcPcscdd/Dee++xbds21qxZg8lkIiMjg7q6Omw2WyQ/oojIsNxuFydOFNDQsJX+s78mIAT4SUkp41/+pYb/+3/hz//cyCpFRGTUWp2Fa0cGsH//fiwWCytXrhxy3el0snr1avx+P9nZ2de8bvLkybz22mvD/s67776b48ePj6jW3Nxcxo0bh81mGwzLRUVFFBUVjej3iUjisVqt+Hw1uFxu6uufIxhMJiWll8zMdDZsqOGJJ6zk5sILL8BjjxldrYhIbPJ6vYOnAJqbm2lubqanpyei72FIn98zZ87gcDiu2d2dNWsWAGfPnh02/I6Wo0ePMnfu3Ki9n4iMTVarlaqq4ce+HToE69fD/fdHuSgRkTgy3ObjqVOnmDdvXsTew5Dw29bWxr333nvN9YkTJw6ui4iMJePGwd//vdFViIhITLc6ExERERGJJEPC76RJk4bd3W1vbx9cFxERERGJNEPC7+zZszl37hyhUGjI9dOnTwOQkZFhRFkiIoZqbW3F6Sxl5sw80tKWMXNmHk5nKa2trUaXJiIyZhgSfvPz8+nq6qK2tnbI9erqamw2G1lZ4UaEioiMTa+80oLdXkB19aMEAi/x1lt1BAIvUl39KNnZBQrAIiIREvEH3g4dOkR3dzednZ1Af+eGyyE3Ly8Ps9nMkiVLyM3NZe3atXR0dGC32/F6vRw5cgSPxxN2MIaIyFi1dWs5nZ3PAws+cNUELKChYSsulztsJwkREbl5EQ+/TzzxBE1NTUD/dLe9e/eyd+9ekpKSaGxsZOrUqQDs27ePsrIynnnmGdrb23E4HOzevZtVq1ZFuiQRkZh34UIAcIdZzaK+/rloliMiMmZFPPw2Njbe1H3jx4+noqKCioqKSJcgIhJ3gsFkINzfepkG1kVE5MNSqzMRkRiQktILhJuMGRpYFxGRD0vhV0QkBmRmpgP+MKv+gXUREfmwFH5FRGKA2+3Cbt8I+IDLbSBDgA+7vQy322VccSIiY4gh441FRGQoq9WKz1eDy+Wmvv45gsFkUlJ6ycxMx+2uwWq1Gl2iiMiYoPArIhIjrFar2pmJiIwyHXsQERERkYSh8CsiEqf++Ed45x2jqxARiS8KvyIicWrdOli4EJqbja5ERCR+KPyKiMSpr34VLlyA+fPhyJFWnM5SZs7MIy1tGTNn5uF0ltLa2mp0mSIiMUUPvImIxKk5c+C11yAvr4UlSwrp63ue/hHJSUCIQKCeEycK8PnULUJE5DLt/IqIxLEpU2DmzPKB4LuAKyOSTcACGhq24nK5jStQRCTGKPyKiMS5kycDQFaY1Szq6wPRLEdEJKYp/IqIxLlgMJkrO75XMw2si4gIKPyKiMS9lJReoC/MamhgXUREQOFXRCTuZWamA/4wq/6BdRERgTgNv//0T//Efffdh8ViIT09ncbGRgAKCwuZMmUKqamppKWlUVlZaXClIiKjz+12YbdvBHxAaOBqCPBht5fhdruMK05EJMbEXauzF198kW3btvHiiy9y//3389vf/paPfvSjAGzatIldu3aRkpLCa6+9xsKFC3nwwQeZPn26wVWLiIweq9WKz1eDy+Wmvv45gsFkUlJ6ycxMx+1WmzMRkQ+Ku/D7zW9+k+9+97vcf//9ANxzzz2Daw6HY/DPycnJpKamYrFYol6jiEi0Wa1WqqrKb3ifzwczZsCkSVEoSkQkBsXVsYfe3l5+8Ytf8Oqrr/KJT3yCe+65h2effXbIPcXFxZjNZnJycti+fTuTJ082qFoRkdgSDMJjj8GCBfDmm0ZXIyJijLgKvxcuXCAYDOLz+Th79iw//elP8Xg87Ny5c/Aej8dDd3c3Xq8Xp9PJ+fPnDaxYRCR2pKTAyy/DuHH9AfiVV4yuSEQk+mI6/Ho8HiwWCxaLhby8PG6//XYA1q9fT2pqKp/85Cd5/PHHOXz48JDXmUwmli9fTk5ODnV1dUaULiISk+65p//oQ2YmLFkC27cbXZGISHRFNPx2dXXhcrlYvHgxVqsVk8nEli1bwt67bt06bDYbZrOZOXPmUFNTM+Se4uJiOjs76ezs5MCBA9x555187GMfu+Z3JSUN39w9GAwyYcKED//BRETGkDvugAMH4PHH+3++/nX4j/9oxeksZebMPNLSljFzZh5OZymtra1GlysiElERDb8XL16ksrKSnp4e8vPzgfDBdMWKFezcuZPNmzdz+PBh5s+fT1FREV6v97rvUVJSgtvtpquri7fffpvKykry8vK4cOECtbW1dHd3EwwG2bNnD36/n9zc3Eh+RBGRMSElBb7/ffje96CiogW7vYDq6kcJBF7irbfqCARepLr6UbKzCxSARWRMiWj4nTZtGpcuXeKnP/0p3/rWt8Led/DgQY4dO8YPf/hD/vqv/5pFixaxfft2cnNzKS0tJRQKhX3tpk2b+NM//VM+/vGPk52dTXFxMatXrwZg27Zt2Gw27rrrLr7//e9TV1eHzWaL5EcUERlT/vt/h89+tpx33nkeWMCVMckmYAENDVtxudzGFSgiEmGj1uqsry/cqE3Yv38/FouFlStXDrnudDpZvXo1fr+f7OzsYV87btw4tm/fzvarDqrdfffdHD9+fES15ubmMm7cOGw222BYLioqoqioaES/T0Qknvz7vweAcAE3i/r656JZjogkMK/XO3gKoLm5mebmZnp6eiL6Hob0+T1z5gwOhwOTaejG86xZswA4e/Zs2PA7Go4ePcrcuXOj9n4iIrEkGEzmyo7v1UwD6yIio2+4zcdTp04xb968iL2HId0e2tramDhx4jXXL19ra2uLdkkiIgkrJaUXCPe3daGBdRGRsSGmW52JiMjoy8xMB/xhVv0D6yIiY4Mh4XfSpEnD7u62t7cProuISHS43S7s9o2AD7j8wHEI8HHPPWXccYcL/YWciIwVhoTf2bNnc+7cuWu6Opw+fRqAjIwMI8oSEUlIVqsVn6+GkpJ9pKcvZcaMZaSnL6WkZB8eTw0ej5XsbPjNb4yuVETkwzPkgbf8/HwqKyupra1l1apVg9erq6ux2WxkZWUZUZaISMKyWq1UVZUPu/bzn8PnP98/EnnfPli0KMrFiYhEUMTD76FDh+ju7qazsxPo79xQW1sLQF5eHmazmSVLlpCbm8vatWvp6OjAbrfj9Xo5cuQIHo8n7GAMERGJPrsdXn0VvvAFyM2Fykr48peNrkpEZGQiHn6feOIJmpqagP7pbnv37mXv3r0kJSXR2NjI1KlTAdi3bx9lZWU888wztLe343A42L1795CdYBERiQ0f/SgcPgxr10JJSf8RiGefBZMemxaROBPx8NvY2HhT940fP56KigoqKioiXYKIiIyCceP6d33T0mD9enjnHfj7vze6KhGRW2PImV8REYlPSUlQWgr33gv33dd/rbW1FZfLTX19gGAwmZSUXjIz03G7XVitVmMLFhG5isKviIjcsvz8/n+2tLTwwAOFNDQ8T/+I5CQgRCBQz4kTBfh8NQrAIhJTdFpLRERGbP368oHgu4ArI5JNwAIaGrbicrmNK05EZBgKvyIiMmL19QEgXHvKrIF1EZHYofArIiIjFgwmc2XH92qmgXURkdih8CsiIiOWktIL9IVZDQ2si4jEDoVfEREZsczMdMAfZtXPrFnp0SxHROSGFH5FRGTE3G4XdvtGwAeEBq6GAB/JyWX827+5COjYr4jEEIVfEREZMavVis9XQ0nJPtLTlzJjxjLS05dSUrKP116r4aMftfLAA/DKK0ZXKiLST31+RUTkQ7FarVRVlQ+79rOfwapVsGQJ/OM/wl/9VZSLExG5inZ+RURk1KSmwksvwVe+0v/zL/9idEUikugUfkVEZFSlpMAPfwgHDsBf/IXR1YhIotOxBxERGXVJSfC5zxldhYiIdn5FREREJIFo51dERAzV2tqKy+Wmvj5AMJhMSkovmZnpuN0urFar0eWJyBgTdzu/EyZMwGKxDP4kJyfz3e9+F4DCwkKmTJlCamoqaWlpVFZWGlytiIhcT0tLC9nZBVRXP0og8BJvvVVHIPAi1dWPkp1dQGtrq9ElisgYE3fht6uri87OTjo7O3nrrbcwmUysWLECgE2bNvH222/T0dHBrl27ePLJJ2lsbDS4YhERCWf9+nIaGp4HFgBJA1dNwAIaGrbicrmNK05ExqS4C78f5PF4eOCBB/jkJz8JgMPhICWl/yRHcnIyqampWCwWI0sUEZHrqK8PAFlhVrMG1kVEIieuw++Pf/xjHnvssSHXiouLMZvN5OTksH37diZPnmxQdSIiciPBYDJXdnyvZhpYFxGJnLgNv7/61a/4zW9+w8qVK4dc93g8dHd34/V6cTqdnD9/3qAKRUTkRlJSeoG+MKuhgXURkciJ6fDr8XgGH2zLy8sbsvbjH/+YZcuWkZqaes3rTCYTy5cvJycnh7q6umiVKyIitygzMx3wh1n1c8896dEsR0QSQETDb1dXFy6Xi8WLF2O1WjGZTGzZsiXsvevWrcNms2E2m5kzZw41NTVD7ikuLh58uO3AgQOD10OhEF6vly996UvXrScYDDJhwoQP/8FERGRUuN0u7PaNgA8IDVwNAT5uu62Ms2ddBIPG1SciY09Ew+/FixeprKykp6eH/Px8AJKShj/LtWLFCnbu3MnmzZs5fPgw8+fPp6ioCK/Xe8P3eeWVV+jp6eHhhx8evHbhwgVqa2vp7u4mGAyyZ88e/H4/ubm5kflwIiIScVarFZ+vhpKSfaSnL2XGjGWkpy+lpGQfb71Vw6FDVlLUkV5EIiii/5Uybdo0Ll26BEBbWxs7duwY9r6DBw9y7NgxvF4vBQUFACxatIimpiZKS0spKCjAZAqfy3ft2kVRUdE192zbto01a9ZgMpnIyMigrq4Om80WoU8nIiKjwWq1UlVVbnQZIpIgRu1/T/f1hXuAAfbv34/FYrnmYTWn08nq1avx+/1kZ2eHff2PfvSja67dfffdHD9+fES15ubmMm7cOGw222BYLioqoqioaES/T0RERERundfrHTwF0NzcTHNzMz09PRF9D0P+MunMmTM4HI5rdm5nzZoFwNmzZ68bfiPt6NGjzJ07N2rvJyIiIiLXGm7z8dSpU8ybNy9i72FIt4e2tjYmTpx4zfXL19ra2qJdkoiIxJk//tHoCkQkHukxAhERiTtvvAHLlsH/+l+t7N/vpr4+QDCYTEpKL5mZ6bjdLqxWq9FlikgMMiT8Tpo0adjd3fb29sF1ERGRcKZPh2nTWnjkkULgecBN/6S4EIFAPSdOFODz1SgAi8g1DDn2MHv2bM6dO0coFBpy/fTp0wBkZGQYUZaIiMSJO++E6dPL6Q++C7gyItkELKChYSsul9uw+kQkdhkSfvPz8+nq6qK2tnbI9erqamw2G1lZWUaUJSIiceTkyQAQ7v9fZFFfH4hmOSISJyJ+7OHQoUN0d3fT2dkJ9HduuBxy8/LyMJvNLFmyhNzcXNauXUtHRwd2ux2v18uRI0fweDxhB2OIiIhcFgwmc2XH92qmgXURkaEiHn6feOIJmpqagP7pbnv37mXv3r0kJSXR2NjI1KlTAdi3bx9lZWU888wztLe343A42L17N6tWrYp0SSIiMgalpPQCfQwfgEMD6yIiQ0U8/DY2Nt7UfePHj6eiooKKiopIlyAiIgkgMzOdQMBP/5nfq/nJzEyPdkkiEgcMOfMrIiLyYbndLuz2jYAPuPwAdQjwYbeX4Xa7jCtORGKW+vyKiEhcslqt+Hw1uFxu6uufu6rPr9qcicjwFH5FRCRuWa1WqqrKjS5DROKIjj2IiEhCeP99CAaNrkJEjKbwKyIiCeHxx+GRR6Cry+hKRMRICr8iIpIQiorgxAlYtAj+4z+MrkZEjKLwKyIiCWHx4v7w+x//AQsWwLlzRlckIkZQ+BURkYTxqU/Bz38OFgs88ADU1bXidJYyc2YeaWnLmDkzD6ezlNbWVqNLFZFRom4PIiKSUD7xCfjZz2Dp0hYeeaQQeB5w0z8pLkQgUM+JEwX4fGqXJjIWaedXREQSzh13wLRp5fQH3wVcGZFsAhbQ0LAVl8ttWH0iMnoUfkVEJCGdPBkAssKsZlFfH4hmOSISJQq/IiKSkILBZK7s+F7NNLAuImONwq+IiCSklJReoC/MamhgXUTGGoVfERFJSJmZ6YA/zKp/YF1ExhqFXxERSUhutwu7fSPgA0IDV0OAD7u9DLfbxbvvGlefiIyOuAu/Z8+eZeHChdx5553Y7XZeeOGFwbXCwkKmTJlCamoqaWlpVFZWGlipiIjEMqvVis9XQ0nJPtLTlzJjxjLS05dSUrIPn6+Gixet2O3w0ktGVyoikRR3fX6//OUvs3z5co4fP84bb7zBokWLyMnJ4f7772fTpk3s2rWLlJQUXnvtNRYuXMiDDz7I9OnTjS5bRERikNVqpaqqfNi18eMhKwseeQR++EP46lejXJyIjIq42/k9d+4cRUVFAMyZMweHw8Gbb74JgMPhICWlP88nJyeTmpqKxWIxrFYREYlft98OtbWwdi08/jh84xvQF+75OBGJG3EXfhcvXsyPf/xjgsEgfr+f3//+92RnZw+uFxcXYzabycnJYfv27UyePNnAakVEJJ4lJ8P3vgff/jY89xw4ndDTY3RVIvJhxF34dbvdVFVVYTab+cxnPsOzzz7LXXfdNbju8Xjo7u7G6/XidDo5f/68gdWKiEi8S0oClws8Hvjnf4a8POjoMLoqERmpmA6/Ho8Hi8WCxWIhLy+Pd955h89+9rN85zvf4f333+ett96ioqKCl656GsFkMrF8+XJycnKoq6szqHoRERlLVq+Gw4fB74cf/ABaW1txOkuZOTOPtLRlzJyZh9NZSmtrq9Glish1RDT8dnV14XK5WLx4MVarFZPJxJYtW8Leu27dOmw2G2azmTlz5lBTUzPknuLiYjo7O+ns7OTAgQOcPXuWYDDIypUrSUpKYvr06SxdupSXX3552PcIBoNMmDAhkh9RREQS2F/+JZw6BV/+cgvZ2QVUVz9KIPASb71VRyDwItXVj5KdXaAALBLDIhp+L168SGVlJT09PeTn5wOQlDT86MgVK1awc+dONm/ezOHDh5k/fz5FRUV4vd6wv99ut9PR0cG+ffvo6+ujqamJn/zkJ8yePZsLFy5QW1tLd3c3wWCQPXv24Pf7yc3NjeRHFBGRBGe3w8aN5TQ0PA8s4MqIZBOwgIaGrbhcbuMKFJHrimirs2nTpnHp0iUA2tra2LFjx7D3HTx4kGPHjuH1eikoKABg0aJFNDU1UVpaSkFBASbTtbl84sSJ7N69m6effpqSkhIsFgurV69mzZo1tLS0sG3bNtasWYPJZCIjI4O6ujpsNlskP6KIiAj19QEgXMDNor7+uWiWIyK3YNT6/PZdpx/M/v37sVgsrFy5csh1p9PJ6tWr8fv9Qzo4fNDnP/95Pv/5z19z/e677+b48eMjqjU3N5dx48Zhs9kGw3JRUdFgSzUREZEPCgaTubLjezXTwLqI3Cqv1zt4CqC5uZnm5mZ6ItxixZAhF2fOnMHhcFyzuztr1iygf4pbuPA7Go4ePcrcuXOj9n4iIhLfUlJ6gT6GD8ChgXURuVXDbT6eOnWKefPmRew9DOn20NbWxsSJE6+5fvlaW1tbtEsSERG5aZmZ6YA/zKqfefPSo1mOiNyCmG51JiIiEovcbhd2+0bAB4QGroYG/r2MxkYXf/iDYeWJyHUYcuxh0qRJw+7utre3D66LiIjEKqvVis9Xg8vlpr7+OYLBZFJSesnMTGf58hq+/GUrCxfCoUPwsY8ZXa2IfJAh4Xf27Nl4vV5CodCQc7+nT58GICMjw4iyREREbprVaqWqqnzYtZ/9DB5+GLKz+wdjOBxRLk5EwjLk2EN+fj5dXV3U1tYOuV5dXY3NZiMrK8uIskRERCIiIwNefRUsFsjJgX/7N6MrEpHLIr7ze+jQIbq7u+ns7AT6OzdcDrl5eXmYzWaWLFlCbm4ua9eupaOjA7vdjtfr5ciRI3g8nrCDMUREROLFJz7RvwP8yCPwhS/Ab38LZrPRVYlIxMPvE088QVNTE9A/3W3v3r3s3buXpKQkGhsbmTp1KgD79u2jrKyMZ555hvb2dhwOB7t372bVqlWRLklERMQQd94JL78Mb72l4CsSKyIefhsbG2/qvvHjx1NRUUFFRUWkSxAREYkZf/InMHu20VWIyGWGPPAmIiKS6FpbWwe6RQSGdItwu11YrVajyxMZsxR+RUREoqylpYUHHiikoeF5wE3/pLgQgUA9J04U4PPVKACLjBINuRAREYmy9evLB4LvAq6MSDYBC2ho2IrL5TauOJExTuFXREQkyurrA0C4tp5ZA+siMhoUfkVERKIsGEzmyo7v1UwD6yIyGhR+RUREoiwlpRfoC7MaGlgXkdGg8CsiIhJlmZnpgD/Mqn9gXURGg8KviIhIlLndLuz2jYAPCA1cDQ38exk2m4u+cBvDIvKhqNWZiIhIlFmtVny+moE+v88N9vmdPz+dyZNr2LrVisMBxcVGVyoy9ij8ioiIGMBqtVJVVT7s2qJFsGRJlAsSSRAKvyIiIjFm6VKjKxAZu3TmV0REREQShsKviIiIiCQMhV8RERERSRg68ysiIhIn+vpg40bIyGjl2DE39fWBwU4RmZnpuN0urFar0WWKxLS42/n91a9+xZ/92Z9xxx13MHPmTP71X/91cK2wsJApU6aQmppKWloalZWVBlYqIiISWb298KtftfDFLxZQXf0ogcBLvPVWHYHAi1RXP0p2dgGtra1GlykS0+Iq/Pb09Pz/9u4/OuryzPv4eyYBG8gECQy/xrLYWGgisCYcApjTgt0GeDYqCQohpFayolvUpbSVQeURAhrbHfSUpVYeIT2wLOmQEEFj+SFisaJmkyq1BWJdG2M40LOQX5IfJZrJzPPHJMGQDCRhku8k83mdM2fJ9/pmvtfce9NzeXPPdZOcnEx6ejqff/4569evJzk5merqagDWr1/P2bNnqa2tZffu3axcuZKysjKDsxYREfGP0FAYOXIT8CwwEzC1RMzATEpLs7DbHYblJ9If9Kvi9+OPP6a2tpYf/vCHmEwmFi9ezOjRo9m/fz8A0dHRhIZ6d3KEhIQQERGBxWIxMmURERG/ev/9EmCGj+gMiotL+jIdkX6nXxW/nk7Oemxubqak5PJf9PT0dMLCwkhISGDbtm2MHDmyL1MUERHpVS5XCJdXfK9kbomLiC/9qvidNGkS4eHhvPDCCzQ1NZGTk0NpaSkNDQ1t9+Tk5NDQ0IDT6SQjI4MzZ84YmLGIiIh/hYY2Ax0Xg7zcLXER8SWgi9+cnBwsFgsWi4WkpCQGDx7M/v372bNnD5o26S0AACAASURBVGPHjuXll19mzpw53HTTTe1+z2w2k5ycTEJCAgUFBQZlLyIi4n/x8TFAkY9oUUtcRHzxa/FbX1+P3W5n7ty5WK1WzGYzGzZs8HnvqlWrsNlshIWFERsbS25ubrt70tPTqauro66ujgMHDgAQGxvLO++8Q2VlJXv37uXTTz8lPj6+02e4XC7Cw8P9+RFFREQM5XDYiYp6EigE3C1X3UAhUVFrcTjsxiUn0g/4tfitrKxk+/btNDU1kZKSAoDJ1Pm+pIULF7Jr1y4yMzM5fPgw06dPJy0tDafTedVnnD59mi+++IK6ujoef/xxxo0bx9y5czl//jz5+fk0NDTgcrnIy8ujqKiIxMREf35EERERQ1mtVgoLc1m2bB8xMXcxceLdxMTcxbJl+ygszFWfX5Fr8OshFxMmTKCmpgaAqqoqsrOzO73v4MGDHD16FKfTSWpqKgCzZ8+mvLyc1atXk5qaitnceV2+c+dOsrOzcbvdzJ07l1deeaUttmXLFpYvX47ZbGby5MkUFBRgs9n8+RFFREQMZ7Va2bFjk9FpiPRLvXbCW2edGVrt378fi8XCokWL2l3PyMhg6dKlFBUVMWvWrE5/d9OmTWza1PEv/OjRo3n77bd7lGtiYiKDBg3CZrO1FctpaWmkpaX16P1EREREpPucTmfbLoBz585x7tw5mpqa/PoMQ443PnXqFNHR0R1Wd6dMmQJ4tzb4Kn57wxtvvEFcXFyfPU9ERKSvfPopvP46rFhhdCYi19bZ4uOJEyeYNm2a355hSLeHqqoqIiMjO1xvvVZVVdXXKYmIiAxIr7wCDz8Ma9bAVf5RViRoGLLyKyIiIn3jJz8Bsxl+/GP43/+F7GwYNMjorESMY0jxO2LEiE5Xd6urq9viIiIi4h+rVsHo0XD//VBRAS++WMGGDQ6Ki0twuUIIDW0mPj4Gh8OubhEy4BlS/E6dOhWn04nb7W637/fkyZMATJ482Yi0REREBqy0NBg5EpKTL/Ctby3hiy+eBRx4j0p2U1JSzPHjqWqXJgOeIXt+U1JSqK+vJz8/v931nTt3YrPZmDFjhhFpiYiIDGiJifDd725qKXxn4i18wVsOzKS0NAu73WFcgiJ9wO8rv4cOHaKhoYG6ujrA27mhtchNSkoiLCyM+fPnk5iYyIoVK6itrSUqKgqn08mRI0fIycnxeTCGiIiIXJ9PPy3Bu+LbmRkUFz/Tl+mI9Dm/F78PP/ww5eXlgPd0t71797J3715MJhNlZWWMHz8egH379rF27VrWrVtHdXU10dHR7Nmzh8WLF/s7JREREWnhcoVwecX3SuaWuMjA5ffit6ysrEv3DR06lM2bN7N582Z/pyAiIiI+hIY2Ax46L4DdLXGRgcuQPb8iIiJijPj4GKDIR7SoJS4ycKn4FRERCSIOh52oqCeBQsDdctUNFBIVtRaHw25cciJ9QIdciIiIBBGr1UphYS52u4Pi4meu6PPrbXP2+uvezhBmLZHJAKTiV0REJMhYrVZ27NjUaezkSfg//we+/3349a91GpwMPPpvOhEREWkzZQr85jewZw8kJ0NDg9EZifiXil8RERFpZ8kSOHAAfv977/aH6mqjMxLxHxW/IiIi0kFiIhw7Bp98At/+Npw9a3RGIv6h4ldEREQ6NX06vPOOd+vD7bfDX/5idEYi10/Fr4iIiPg0aRK8+y5ERnqL34qKCjIyVnPrrUlMmnQ3t96aREbGaioqKoxOVaRL1O1BRERErspmg/ffh+rqC8yatYTS0mcBB95T4tyUlBRz/HgqhYXeVmkigUwrvyIiInJNoaGwZs2mlsJ3JpePRzYDMyktzcJudxiXoEgXqfgVERGRLikuLgFm+IjOaImLBDYVvyIiItIlLlcIl1d8r2RuiYsENhW/IiIi0iWhoc2Ax0fU3RIXCWwqfkVERKRL4uNjgCIf0SKmTYvpy3REeiQgi9+tW7cSFxfH4MGD2bBhQ7tYRUUFSUlJhIeHM2nSJI4ePdouvmTJEsaMGUNERASTJk1i+/btfZm6iIjIgOVw2ImKehIoBNwtV90tP6/ls8/s1NUZlp5IlwRk8Ttu3Dg2btxIcnIyJlP7vUWPPPII48aNo7Kykueee47FixdT/ZVzF9evX8/Zs2epra1l9+7drFy5krKysr7+CCIiIgOO1WqlsDCXZcv2ERNzFxMn3k1MzF0sW7aPV17J5U9/svLd74Ja/kogC8g+vwsWLADg1VdfxeO5vLeovr6eV199lbKyMr72ta9x11138Y//+I+8+uqrZGRkABAdHd12f0hICBEREVgslr79ACIiIgOU1Wplx45Nncbeegvmz/ceh3zkCIwf37e5iXRFQK78+vLJJ58QHh7OuHHj2q5NmTKF06dPt7svPT2dsLAwEhIS2LZtGyNHjuzrVEVERIJObKz3NLgvv/Qeh1yizmcSgPpV8VtfX09ERES7axEREdTX17e7lpOTQ0NDA06nk4yMDM6cOdOXaYqIiAStW265fBzyT39qdDYiHRle/Obk5GCxWLBYLCQlJV313vDwcGpra9tdu3jxYqfbGsxmM8nJySQkJFBQUODXnEVERMS3sWPh7bfhv/7L6ExEOup28VtfX4/dbmfu3LlYrVbMZnOHjgxfvXfVqlXYbDbCwsKIjY0lNze33T3p6enU1dVRV1fHgQMHOrzHV7/w9s1vfpP6+nr+9re/tV07efIkt956q898XS4X4eHh3f2YIiIich1uvBG061ACUbeL38rKSrZv305TUxMpKSkAHToytFq4cCG7du0iMzOTw4cPM336dNLS0nA6nVd9RnNzM42NjbhcLpqammhsbMTtdhMeHs6CBQtYv349jY2NvPbaa/z5z3/m7rvvBuD8+fPk5+fT0NCAy+UiLy+PoqIiEhMTu/sxRURERGQA6na3hwkTJlBTUwNAVVUV2dnZnd538OBBjh49itPpJDU1FYDZs2dTXl7O6tWrSU1NxWzuvPZ++umn2bhxY9vPWVlZ7Ny5kx/84Ae8+OKL3H///YwYMYKbbrqJvLw8IiMj2+7dsmULy5cvx2w2M3nyZAoKCrDZbN39mCIiItILKioqsNsdFBeX4HKFEBraTHx8DA6HHavVanR6EgRMnq/2EuumyspKRo0aRWZmJuvWrWsXe/DBB8nLy6OmpqZdkbtnzx6WLl3Ku+++y6xZs3qeuR+cOHGCadOmERkZyaBBg7DZbG2FclpaGmlpaYbmJyIiMpBcuHCB229fQmnps8AMwIT3kIxioqKepLAwVwVwkHM6nW07BM6dO8e5c+doamqiurqaDz74gLi4uOt+Rq/1+T116hTR0dEdVnenTJkCwOnTpw0vflu98cYbfhlMERER8W3Nmk0the/Mr1w1AzMpLc3Cbnf47CEswaGzxcfWxUp/6bVuD1VVVe22I7RqvVZVVdVbjxYREZEAVFxcgnfFtzMzWuIivcvwVmciIiISHFyuELxbHTpjbomL9K5eK35HjBjR6epudXV1W1xERESCR2hoM+Drq0bulrhI7+q14nfq1Kl89NFHuN3udtdPnjwJwOTJk3vr0SIiIhKA4uNjgCIf0SK+9a2YvkxHglSvFb8pKSnU19eTn5/f7vrOnTux2WzMmOFrz4+IiIgMRA6HnaioJ4FCvF0eaPm/hQwatJZ33rFz8aJx+Ulw6FG3h0OHDtHQ0EBdXR3g7dzQWuQmJSURFhbG/PnzSUxMZMWKFdTW1hIVFYXT6eTIkSPk5OT4PBhDREREBiar1UphYW5Ln99n2vX5tdtz+eMfrQwbZnSWMtD1qM/vzTffTHl5ufcNTCZa38JkMlFWVsb48eMBaGhoYO3ateTl5VFdXU10dDRPPPEEixcv9uNH6LnW1hn+6hsnIiIiIv7l73qtRyu/ZWVlXbpv6NChbN68mc2bN/fkMSIiIiIifqVWZyIiIiISNFT8ioiIiEjQUPErIiIiAa2uDp57Dq7onirSIz3a8ysiIiLSV37/e7Db4b33KggPd/DBByXtOkU4HHasVqvRaUo/oeJXREREAtqdd8L27RdYvnwJ8CzgwHtMspuSkmKOH0+lsDBXBbB0ibY9iIiISMB7551NeAvfmXgLX/CWMTMpLc3CbncYlpv0Lyp+RUREJOAVF5cAvk6HndESF7k2Fb8iIiIS8FyuEC6v+F7J3BIXuTYVvyIiIhLwQkObAV+H0rpb4iLXpuJXREREAl58fAxQ5CNa1BIXuTYVvyIiIhLwHA47UVFPAoVAa8NfN1BIVNRaHA67cclJv6JWZyIiIhLwrFYrhYW52O0OioufuaLPr9qcSdep+BUREZF+wWq1smPHJqPTkH5O2x5ERERkwPnrX43OQAKVil8REREZUEpL4dZb4f/+X/D4ahAhQSsgi9+tW7cSFxfH4MGD2bBhQ5djAEuWLGHMmDFEREQwadIktm/f3ldpi4iISACIioKsLO/r0UfB7b7270jwCMjid9y4cWzcuJHk5GRMJlOXYwDr16/n7Nmz1NbWsnv3blauXElZWVlfpS4iIiIB4LHHIDsb/t//g/vug6YmozOSQBGQX3hbsGABAK+++iqeK/694moxgOjo6LY/h4SEEBERgcVi6cVsRUREJBA98AAMGwZLl8LFi/CrX1WQmemguLjkim4RdnWLCCIBWfxer/T0dPbt2wfAnj17GDlypMEZiYiIiBHuvddbAC9YcIFvfWsJjY3PAg68RyW7KSkp5vjxVAoL1S4tWATktofrlZOTQ0NDA06nk4yMDM6cOWN0SiIiImKQxES4445NLYXvTLyFL3jLoJmUlmZhtzuMS1D6lOHFb05ODhaLBYvFQlJSkt/e12w2k5ycTEJCAgUFBX57XxEREel/PvusBJjhIzqD4uKSvkxHDNTt4re+vh673c7cuXOxWq2YzeZOuy603rtq1SpsNhthYWHExsaSm5vb7p709HTq6uqoq6vjwIEDHd6jsy+1dSXWyuVyER4efs37REREZOByuUK4vOJ7JXNLXIJBt4vfyspKtm/fTlNTEykpKYDvInThwoXs2rWLzMxMDh8+zPTp00lLS8PpdF71Gc3NzTQ2NuJyuWhqaqKxsRF3S5+Sq8XOnz9Pfn4+DQ0NuFwu8vLyKCoqIjExsbsfU0RERAaQ0NBmwFfTX3dLXIJBt4vfCRMmUFNTw7Fjx/jZz37m876DBw9y9OhRtm7dyoMPPsjs2bPZtm0biYmJrF69uq1g7czTTz/NkCFD+M///E+ysrIYMmQIu3fvvmYMYMuWLdhsNkaNGsULL7xAQUEBNputux9TREREBpD4+BigyEe0qCUuwcDk6axfWBdVVlYyatQoMjMzWbduXbvYgw8+SF5eHjU1NZjNl2vsPXv2sHTpUt59911mzZrV88z94MSJE0ybNo3IyEgGDRqEzWZrK5TT0tJIS0szND8RERHxj4qKCmbNSqW0NAvv3l8z4AaKiIpay3vv5TJqlLo9GM3pdLbtEDh37hznzp2jqamJ6upqPvjgA+Li4q77Gb3W6uzUqVNER0e3K3wBpkyZAsDp06cNL35bvfHGG34ZTBEREQlMVquVwsJc7HYHxcXPtOvzm5WVyyOPWFm0CBYvNjrT4NbZ4mPrYqW/9FrxW1VVxS233NLhemRkZFtcREREpK9YrVZ27NjU4XpTE9xwAyxZ4j0M48EHDUhO+syAPORCREREpKsGDYJdu7yHYTz0kLcAfuwxo7OS3tJrxe+IESM6Xd2trq5ui4uIiIgEArMZXngBbrwRVq+Gmhp45hnoQldV6Wd6rfidOnUqTqcTt9vdbt/vyZMnAZg8eXJvPVpERESk20wmyMqC4cO9BfDnn8Mvf+ktjGXg6LX/d6akpFBfX09+fn676zt37sRmszFjhq9TVkRERESM89hjsG0bbN0KP/2p0dmIv/Vo5ffQoUM0NDRQV1cHeDs3tBa5SUlJhIWFMX/+fBITE1mxYgW1tbVERUXhdDo5cuQIOTk5XTqdTURERMQIDz7oXQH+1reMzkT8rUfF78MPP0x5eTngPd1t79697N27F5PJRFlZGePHjwdg3759rF27lnXr1lFdXU10dDR79uxhsfqIiIiISIC7997Lf66oqGhpk1bSrk2aw2HHalV/4P6kR8VvWVlZl+4bOnQomzdvZvPmzT15jIiIiIjhLly4wO23L6G09FnAAZgANyUlxRw/nkphYa4K4H5EW7hFRERErmLNmk0the9MvIUveEuomZSWZmG3O4xLTrpNxa+IiIjIVRQXl+A9ErkzM1ri0l+o+BURERG5CpcrhMsrvlcyt8Slv1DxKyIiInIVoaHNgMdH1N0Sl/5Cxa+IiIjIVcTHxwBFPqJFjB8f05fpyHVS8SsiIiJyFQ6HnaioJ4FCwN1y1Q0U8rWvreXYMTsHDxqXn3SPil8RERGRq7BarRQW5rJs2T5iYu5i4sS7iYm5i2XL9vE//5PLvHlWFiyAvDyjM5Wu6FGfXxEREZFgYrVa2bFjU6ex/HzIyIC0NKirgwce6OPkpFu08isiIiJyHQYNgl274KGH4J13wOPru3ESELTyKyIiInKdzGZ48UVwu8HkqyuaBASt/IqIiIj4gckEIWr5G/BU/IqIiIhI0NC2BxEREZFeVlFRgd3uoLi4BJcrhNDQZuLjY3A47FitVqPTCyoqfkVERER60YULF7j99iWUlj4LOPAeleympKSY48dTKSzMVQHch7TtQURERKQXrVmzqaXwnYm38AVvCTaT0tIs7HaHcckFoYAsfrdu3UpcXByDBw9mw4YNbde//PJLMjIyGD9+PMOGDWPWrFkUFha2+90lS5YwZswYIiIimDRpEtu3b+/r9EVERETaFBeXADN8RGe0xKWvBGTxO27cODZu3EhycjKmr/QLcblcfOMb3+C9997j4sWLrFixgrvvvptLly613bN+/XrOnj1LbW0tu3fvZuXKlZSVlRnxMURERERwuUK4vOJ7JXNLXPpKQBa/CxYs4M4772TYsGF4vtIpesiQITz11FPcdNNNAPzgBz/A7Xbz17/+te2e6OhoQkO9W5lDQkKIiIjAYrH07QcQERERaREa2gz4OvnC3RKXvhKQxW9X/eUvf+HSpUtERUW1u56enk5YWBgJCQls27aNkSNHGpShiIiIBLv4+BigyEe0qCUufaXfFr9///vfue+++3jqqacYMmRIu1hOTg4NDQ04nU4yMjI4c+aMQVmKiIhIsHM47ERFPQkUAu6Wq26gkNDQtWRl2Y1LLggZXvzm5ORgsViwWCwkJSV16XeamppYtGgRkydP5oknnuj0HrPZTHJyMgkJCRQUFPgzZREREZEus1qtFBbmsmzZPmJi7mLixLuJibmLJUv28eqruYwbpzZnfanbxW99fT12u525c+ditVoxm83tOjJcee+qVauw2WyEhYURGxtLbm5uu3vS09Opq6ujrq6OAwcOdHgP0xUHZLvdbu677z4GDx7Mr3/962vm63K5CA8P78YnFBEREfEvq9XKjh2bOH36AB9/XMDp0wdwOjfxz/+swrevdbv4raysZPv27TQ1NZGSkgJ0LFBbLVy4kF27dpGZmcnhw4eZPn06aWlpOJ3Oqz6jubmZxsZGXC4XTU1NNDY24nZ7/5ngX//1Xzl//jx79uzBbG6f/vnz58nPz6ehoQGXy0VeXh5FRUUkJiZ292OKiIiIyADU7RPeJkyYQE1NDQBVVVVkZ2d3et/Bgwc5evQoTqeT1NRUAGbPnk15eTmrV68mNTW1Q/Ha6umnn2bjxo1tP2dlZbFz507mzJnDr3/9a8LCwtp9ie3w4cMkJCQAsGXLFpYvX47ZbGby5MkUFBRgs9m6+zFFREREZAC6ruONv9qG7Er79+/HYrGwaNGidtczMjJYunQpRUVFzJo1q9PfzczMJDMzs9NY6wpwZ0aPHs3bb7997cSvkJiYyKBBg7DZbG2FclpaGmlpad1+LxERERHpGafT2bZD4Ny5c5w7d46mpia/PuO6it+rOXXqFNHR0R1Wd6dMmQLA6dOnfRa/fe2NN94gLi7O6DREREREAPjkE3jzTfjhD43OpG91tvh44sQJpk2b5rdn9Fq3h6qqKiIjIztcb71WVVXVW48WERER6ddeeQVWrICNG+Eq/9AuPdBrK78iIiIi0jOrV4PLBU8+Cf/7vxX8/e8O/vCHElyuEEJDm4mPj8HhsGO1qltEd/Va8TtixIhOV3erq6vb4iIiIiLSOe9RBhd48sklwLOAAzABbkpKijl+PJXCwlwVwN3Ua9sepk6dykcffdThC2onT54EYPLkyb31aBEREZEB4X/+ZxPewncm3sIXvOXbTEpLs7DbHYbl1l/1WvGbkpJCfX09+fn57a7v3LkTm83GjBkzeuvRIiIiIgNCcXEJ4KtmmtESl+7o0baHQ4cO0dDQQF1dHeDt3NBa5CYlJREWFsb8+fNJTExkxYoV1NbWEhUVhdPp5MiRI+Tk5Pg8GENEREREvFyuEC6v+F7J3BKX7uhR8fvwww9TXl4OeE9327t3L3v37sVkMlFWVsb48eMB2LdvH2vXrmXdunVUV1cTHR3Nnj17WLx4sf8+gYiIiMgAFRraDHjovAB2t8SlO3pU/JaVlXXpvqFDh7J582Y2b97ck8eIiIiIBLX4+BhKSorw7vm9UhHx8TF9nVK/12t7fkVERETk+jgcdqKingQKgdYmAm6gkKiotTgcduOS66fU51dEREQkQFmtVgoLc7HbHRQXP3NFn1+1OesJFb8iIiIiAcxqtbJjx6Zr3lddDZ0critX0LYHERERkX4uOxtiYuDUKaMzCXwqfkVERET6ubvvhrFjYfZseP99o7MJbCp+RURERPq5UaPgd7+DiRPhu9+Fd94xOqPApeJXREREZAAYPhyOHIFp02DuXHjjDaMzCkwqfkVEREQGCIsFDh6EOXPgzjuhoMDojAKPuj2IiIiIDCBhYfDKK7B0KaSmwh/+UMHzzzsoLi65olWaPShbpan4FRERERlgBg+GPXvg8OELJCcvobT0WcCB95hkNyUlxRw/nkphYfD1Cta2BxEREZEBKDQUXn55U0vhOxNv4Qve8m8mpaVZ2O0O4xI0iIpfERERkQGquLgEmOEjOqMlHlxU/IqIiIgMUC5XCJdXfK9kbokHFxW/IiIiIgNUaGgz4PERdbfEg4uKXxEREZEBKj4+BijyES1i+vSYvkwnIARk8bt161bi4uIYPHgwGzZsaBdbsmQJY8aMISIigkmTJrF9+/ZuxUVERESChcNhJyrqSaAQcLdcdbf8vJaGBjsul2HpGSIgi99x48axceNGkpOTMZna71NZv349Z8+epba2lt27d7Ny5UrKysq6HBcREREJFlarlcLCXJYt20dMzF1MnHg3MTF3sWzZPrZty+WVV6wsWQJffml0pn0nIPv8LliwAIBXX30Vj6f9PpXo6Oi2P4eEhBAREYHFYulyXERERCSYWK1WduzY1GlszBi4915ISYH8fO8BGQNdQK78Xkt6ejphYWEkJCSwbds2Ro4c2a24iIiIiMBdd8FvfwvHjkFSEtTXG51R7+uXxW9OTg4NDQ04nU4yMjI4c+ZMt+IiIiIi4pWYCK+/Du+/D/PmwcWLRmfUuwwvfnNycrBYLFgsFpKSkrr8e2azmeTkZBISEigoKOh2XERERES8vv1tePNNuOEGaB7g3c+6XfzW19djt9uZO3cuVqsVs9ncoSPDV+9dtWoVNpuNsLAwYmNjyc3NbXdPeno6dXV11NXVceDAgQ7vceUX3q7kcrkIDw/vcVxEREREYPp0bwEcGWl0Jr2r28VvZWUl27dvp6mpiZSUFMB3gbpw4UJ27dpFZmYmhw8fZvr06aSlpeF0Oq/6jObmZhobG3G5XDQ1NdHY2Ijb7eb8+fPk5+fT0NCAy+UiLy+PoqIiEhMTAa4ZFxERERHfrrHmOCB0u9vDhAkTqKmpAaCqqors7OxO7zt48CBHjx7F6XSSmpoKwOzZsykvL2f16tWkpqZiNndeez/99NNs3Lix7eesrCx27tzJvHnz2LJlC8uXL8dsNjN58mQKCgqw2Wxt914rLiIiIiJdU1FRgd3uoLi4BJcrhNDQZuLjY3A47FitVqPT65HranV2ZRuyr9q/fz8Wi4VFixa1u56RkcHSpUspKipi1qxZnf5uZmYmmZmZncbefvttn88cPXr0VeO+JCYmMmjQIGw2W1uhnJaWRlpaWrffS0RERGQguHDhArffvoTS0mcBB2AC3JSUFHP8eCqFhbl+L4CdTmfbDoFz585x7tw5mpqa/PqMXuvze+rUKaKjozus7k6ZMgWA06dP+yx++9obb7xBXFyc0WmIiIiIBIw1aza1FL4zv3LVDMyktDQLu93hs39wT3W2+HjixAmmTZvmt2f0WreHqqoqIjvZMd16raqqqrceLSIiIiLXqbi4BJjhIzqjJd7/GN7qTEREREQCj8sVgnerQ2fMLfH+p9eK3xEjRnS6ultdXd0WFxEREZHAFBraDPj6fpe7Jd7/9FrxO3XqVD766CPcbne76ydPngRg8uTJvfVoEREREblO8fExQJGPaBFjx8b0ZTp+02vFb0pKCvX19eTn57e7vnPnTmw2GzNm+NpDIiIiIiJGczjsREU9CRQCrYuZbqCQoUPX8tZbdq4o8/qFHnV7OHToEA0NDdTV1QHezg2tRW5SUhJhYWHMnz+fxMREVqxYQW1tLVFRUTidTo4cOUJOTs41T24TEREREeNYrVYKC3Nb+vw+067Pb1ZWLo89ZmXJEjh9GiZNMjrbrjN5rtas14ebb76Z8vJy7xuYTG39fk0mE2VlZYwfPx6AhoYG1q5dS15eHtXV1URHR/PEE0+wePFiP36EnmttnfHBBx+o1ZmIiIhINzQ3w+uvwz//c+8+x9/1Wo9WfsvKyrp039ChQ9m8/xQITAAADvBJREFUeTObN2/uyWNEREREJECFhPR+4dsb1OpMRERERIKGil8RERERCRoqfkVEREQkaKj4FRERERG/+vvf4Ze/hCuOewgIPfrCm4iIiIiIL7/7HfzoR/Df/11BaKiD998vadcqzeGwY7VaDclNxa+IiIiI+NWdd8Ivf3mBRx9dAjwLOAAT4KakpJjjx1MpLMw1pADWtgcRERER8bv339+Et/CdibfwBW/pOZPS0izsdochean4FRERERG/Ky4uAWb4iM5oifc9Fb8iIiIi4ncuVwiXV3yvZG6J9z0VvyIiIiLid6GhzYDHR9TdEu97Kn5FRERExO/i42OAIh/RopZ431PxKyIiIiJ+53DYiYp6EigEWhv+uoFCoqLW4nDYDclLrc5ERERExO+sViuFhbnY7Q6Ki5+5os+vMW3OQMWviIiIiPQSq9XKjh2bjE6jHW17EBERERHDNffR998CsvjdunUrcXFxDB48mA0bNnR6T2FhIWazmaysrHbXlyxZwpgxY4iIiGDSpEls3769L1IWERERkR565x2YOhU+/bT3nxWQxe+4cePYuHEjycnJmEwd+8O53W5+/OMfM2vWrA7x9evXc/bsWWpra9m9ezcrV66krKysr1IXERERkW6aMAGamuA734GPP+7dZwVk8btgwQLuvPNOhg0bhsfTsT/cSy+9REJCAhMnTuwQj46OJjTUu5U5JCSEiIgILBZLn+QtIiIiIt13003w+9/DsGEwezacOtV7zwrI4vdqqqqq2LJlC+vXr/d5T3p6OmFhYSQkJLBt2zZGjhzZhxkGB6fTaXQK/Y7GrGc0bt2nMesZjVv3acx6RuPWubFj4a23YMwYmDMH/vjH3nlOvyt+n3jiCX76058SEREB0Om2iJycHBoaGnA6nWRkZHDmzJm+TnPA01/c7tOY9YzGrfs0Zj2jces+jVnPaNx8s1rhd7+Db3wD5syp4M47V3PvvSv9+gzDi9+cnBwsFgsWi4WkpKSr3vvBBx/wxz/+kQceeAAAj8fT6bYIALPZTHJyMgkJCRQUFPg9bxERERHxv8hI+M1vLvDFF6kcOHAPZWX/4df373bxW19fj91uZ+7cuVitVsxms8+ODPX19axatQqbzUZYWBixsbHk5ua2uyc9PZ26ujrq6uo4cOBAh/f46sruO++8Q0lJCaNGjcJqtZKbm8vPfvYzli1b5jNfl8tFeHh4dz+miIiIiBgkK2sTX3zxLDAT6Piv/Nej28VvZWUl27dvp6mpiZSUFKDzrQcACxcuZNeuXWRmZnL48GGmT59OWlraNZf7m5ubaWxsxOVy0dTURGNjI263m+XLl/PJJ5/wpz/9iQ8//JC7776bRx99lF/84hcAnD9/nvz8fBoaGnC5XOTl5VFUVERiYmJ3P6aIiIiIGKS4uASY0Svv3e0T3iZMmEBNTQ3g/fJZdnZ2p/cdPHiQo0eP4nQ6SU1NBWD27NmUl5ezevVqUlNTMZs7r72ffvppNm7c2PZzVlYWO3fu5Ac/+AFDhw5tuz5kyBAiIiIYPnx427UtW7awfPlyzGYzkydPpqCgAJvN1ulzLl26BMBHH33UjREQgIsXL3LixAmj0+hXNGY9o3HrPo1Zz2jcuk9j1jMat2urr68DWr/x5q3TWuu26+a5DhUVFR6TyeTZsGFDh9jy5cs9ERERnubm5nbXnU6nx2Qyed57773rebRf7N692wPopZdeeumll1566RXgr927d/ul/uv2ym9XnTp1iujo6A6ru1OmTAHg9OnTzJo1q7ce3yXz5s1j9+7dTJgwgbCwMENzEREREZGOLl26xGeffca8efP88n69VvxWVVVxyy23dLgeGRnZFjfayJEjSU9PNzoNEREREbmKhIQEv72X4a3ORERERET6Sq8VvyNGjOh0dbe6urotLiIiIiLSl3qt+J06dSofffQRbre73fWTJ08CMHny5N56tIiIiIhIp3qt+E1JSaG+vp78/Px213fu3InNZmPGjN7p3SYiIiIi4kuPit9Dhw6Rn5/Pa6+9Bng7N+Tn55Ofn9/Wg23+/PkkJiayYsUKsrOzOXbsGA899BBHjhzB4XD4PBijL3Tl5Dlp76233sJsNnf6Ki4uNjo9w3Xn5MMTJ07wve99D4vFwvDhw7nnnnsoKyvr44wDQ1fHbdmyZZ3OvZiYGAOyNtabb77J/fffz8SJExk6dCg33XQTycnJnfYM1Vy7rKvjprl22YcffkhSUhL/8A//wJAhQxgxYgS33347OTk5He7VXLusq+OmuXZ12dnZmM1mLBZLh9j1zrcedXt4+OGHKS8vB7ynu+3du5e9e/diMpkoKytj/PjxAOzbt4+1a9eybt06qquriY6OZs+ePSxevLgnj/WbhQsX8v777/Pv//7vTJw4kZycHNLS0nC73aSlpRmaW6D72c9+xh133NHu2q233mpQNoGj9eTD2267jZSUFLKzszv9D7y//OUvzJkzh7i4OPbu3culS5dYt24d3/72t/nwww8ZOXKkAdkbp6vjBhAWFsaxY8c6XAs2L730EhUVFfz4xz/m1ltvpaKigueff56ZM2fy+uuvt/391Fxrr6vjBpprrS5evMj48eNJT0/HZrNRX19PTk4O9913H5999hlr164FNNeu1NVxA801X86dO8djjz3GuHHjqK2tbRfzy3zzS7fgfuTAgQMek8nk2bNnT7vrc+fO9dhstg6HcojXsWPHPCaTyfPyyy8bnUrAq6ys9Hn4y6JFizyjRo3y1NXVtV0rLy/3DB482LNmzZq+TDPgXG3c7r//fo/FYjEgq8Bz/vz5Dtfq6+s9Y8aM8Xzve99ru6a51l5Xx01z7dpmzpzpGT9+fNvPmmtdc+W4aa75duedd3qSk5M9y5Yt84SHh7eL+WO+BV2rs/3792OxWFi0aFG76xkZGfztb3+jqKjIoMz6B4/HY3QKAc/XGLlcLn77299yzz33EB4e3nZ9/Pjx3HHHHezfv7+vUgxI15pbmnteo0aN6nBt6NChREdHc/bsWUBzrTNdGbdWmmtXN2LECEJDvf9wrLnWdV8dt1aaax3t3r2b48eP86tf/arD+PhrvgVd8duVk+fEt0ceeYRBgwYxbNgw5s+fz7vvvmt0Sv1GaWkpjY2NTJ06tUNsypQp/PWvf+XLL780ILP+4dKlS4wdO5bQ0FC+/vWv82//9m/U1NQYnVZAuHjxIidOnGjbgqS51jVXjlsrzbX2PB4PLpeLiooKXnzxRV5//XUee+wxQHPtaq42bq0019o7f/48q1at4uc//znjxo3rEPfXfOu1E94CVX84eS4Q3XjjjaxatYo5c+YwYsQIPvnkEzZt2sScOXM4cOAAc+fONTrFgNc6t1rn2ldFRkbi8Xioqalh9OjRfZ1awLvtttuIjY1ta5H41ltv8Ytf/II333yTP/zhDwwdOtTgDI31yCOPcOnSpba9hJprXXPluIHmWmdWrFjBtm3bAAgJCeG5555jxYoVgOba1Vxt3EBzrTOPPPIIMTEx/PCHP+w07q/5FnTFr/TMbbfdxm233db2c0JCAikpKUyZMoU1a9ao+JVetWrVqnY//9M//ROxsbHce++9ZGdn86Mf/cigzIz31FNP8Zvf/IYXXniB2NhYo9PpN3yNm+ZaR2vXruWhhx7iwoULFBQU8JOf/ITGxkbWrFljdGoB7VrjprnWXn5+Pr/97W/505/+1OvPCrriVyfP+c+wYcNISkripZde4osvvuCGG24wOqWA1jq3WufaV1VXV2MymRg+fHhfp9VvpaSkMHTo0KDep79hwwaysrJ49tlnefjhh9uua65dna9x8yXY59rXv/51vv71rwPeNqbg/Y+HBx54QHPtKnyN27/8y79gtVo7/Z1gnWv19fU8+uijrFy5ktGjR/P5558DtG1huHjxIqGhoX6bb0G351cnz/UOI/s29xdRUVGEhYXx5z//uUPs5MmTfPOb32Tw4MEGZNY/eTyeDn+Pg8mGDRvaXo8//ni7mOaab1cbN1+Cfa5dafr06bhcLj799FPNtW5oHber9aMN1rlWWVnJhQsXeO6554iMjGx77dmzh4aGBoYPH859993HLbfc4pf5FnTFr06e85+amhpee+01YmNj9T9uXRAaGspdd93Fvn37qK+vb7t+5swZjh07xsKFCw3Mrv9pPVRn1qxZRqfS555++mk2bNjAU089xVNPPdUhrrnWuWuNmy/BPNc6c+zYMUJCQoiKitJc64avjpsvwTrXxo4dy7Fjx3jrrbfaXseOHWPevHl87Wtf46233uKZZ54hJCTEL/PN5AnCPhvz5s1rO+QiKioKp9NJdnZ222EX0lF6ejo333wzcXFxREZG8sknn/D8889TVlbGoUOH+O53v2t0ioY7dOgQDQ0N1NXV8cADD7Bo0aK2lnpJSUmEhYXx8ccfM336dOLi4nj88cfbmnN//vnnfPjhh0G57eZa43bhwgW+//3vs3TpUr7xjW/g8Xj4/e9/z3/8x39wyy23UFRUFFRN4Z9//nlWr17N/PnzWb9+fYdWQDNnzgTQXLtCV8atvLxcc+0rHnroIYYNG8b06dMZPXo0lZWV7N27l7y8POx2Oz//+c8BzbUrdWXcNNe6ZtmyZbz88svU1dW1XfPLfOtZ++H+rb6+3vOjH/3IM3bsWM8NN9zgue222zy5ublGpxXQfv7zn3tiY2M9N954oyc0NNQzatQozz333ON5//33jU4tYEyYMMFjMpk8JpPJYzab2/25vLy87b4PPvjA873vfc8zdOhQz7BhwzwLFy70fPrppwZmbqxrjVtNTY1n4cKFnptvvtkzZMgQzw033OCZNGmS5/HHH/fU1tYanX6fmzNnTrtx+urLbDa3u1dz7bKujJvmWns7duzwfOc73/FYrVbPoEGDPMOHD/fccccdnpycnA73aq5d1pVx01zrmmXLlnV6EMj1zregXPkVERERkeAUdHt+RURERCR4qfgVERERkaCh4ldEREREgoaKXxEREREJGip+RURERCRoqPgVERERkaCh4ldEREREgoaKXxEREREJGip+RURERCRoqPgVERERkaCh4ldEREREgsb/BzpNlyA7bXC0AAAAAElFTkSuQmCC",
      "text/plain": [
       "PyPlot.Figure(PyObject <matplotlib.figure.Figure object at 0x7f784acfda10>)"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "1-element Array{Any,1}:\n",
       " PyObject <matplotlib.lines.Line2D object at 0x7f784acc7c50>"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "semilogy(Q, err, \"o--\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "calc_deriv_err (generic function with 3 methods)"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "function calc_deriv_err{T<:Number, QT<:Jacobi.QUADRATURE_TYPE}(Q, ::Type{QT}=Jacobi.GLJ, ::Type{T}=Float64)\n",
    "    \n",
    "    z = Jacobi.qzeros(QT, Q, 0, 0, T)\n",
    "    w = Jacobi.qweights(QT, z, 0, 0)\n",
    "    D = Jacobi.qdiff(QT, z)\n",
    "    \n",
    "    y = fun(z)\n",
    "    yex = dfun(z)\n",
    "    return maxabs(D * y - yex)\n",
    "end\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "Q = 2:40\n",
    "tipo =BigFloat\n",
    "QT = Jacobi.GLJ\n",
    "derr = [calc_deriv_err(q, QT, tipo)  for q in Q];"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAArgAAAIICAYAAACFCp1aAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3X9w0/ed5/GnZJOtIYIUEEmjrEsiGtc/YAYYbAxt2LvELrcGismBMW4S1KFNYXpZkq3FNk4IUOB2RW7jsruXNiZrl8QrG7tmQ/kVoHez8bU+eSe0PUAk9FzXSdwrFjaJfyQpluX7Q8bEIBlsLH9l+fWYyYR8P/ro+5YmwGve/nw/H1Nvb28vIiIiIiIxwmx0ASIiIiIiI0kBV0RERERiigKuiIiIiMQUBVwRERERiSkKuCIiIiISUxRwRURERCSmKOCKiIiISEyJN7qASLt06RJvvvkmM2fOJCEhwehyREREROQ6n3zyCb///e/52te+xvTp02/7/WI+4L755pt84xvfMLoMEREREbmJ119/nYKCgtt+n5gPuDNnzgSCX1hycnLE75ebu5n33isOO56YuJmDB8OPR5Onn36al156yegyxhR9Z8Oj723o9J0Nj763odN3Njz63obm/PnzfOMb3+jPbbcr5gPu1WUJycnJzJs3L+L3u/NOCzAXMIUYDWAyWUaljpEwZcqUMVNrtNB3Njz63oZO39nw6HsbOn1nw6PvbXhGajmpHjIbYenpKYAnzKiHpqYUcnLg3LnRrEpERERk/FDAHWEulxO7/VmgDgj0XQ0AddjtRezb5+Sdd2DOHPjnfzauThEREZFYFfNLFEab1Wqlrq4Sp9NFff1O/P444uN7SE9PweWqxGq18thj8PLL8PDDRlcrIiIiEnsUcCPAarVSWron7Pgdd8Bf/VX4+T6fry8ge68LyE6sVuug976dudfLz88f0utF39lw6XsbOn1nw6Pvbej0nQ2PvjdjmXp7e3uNLiKSTp8+zfz583n77bfHxGLvlpYWFi1aS0PDbiCD4MNqAaAeu/1Z6uoqwwbV25krIiIiYpSRzmtagxtltmzZ0xdQF3JtJwYzsJCGhl04na6IzBURERGJFQq4Uaa+3kuw+xpKBm63l69+FTZuHPrc4LiIiIhIbFPAjTJ+fxyh99AFMJOQEMcDD8DkyUOfGxwXERERiW16yCzKxMf3AL2EOyji3nt7+MlPhjc3OC4iIiIS29TBjTI3OygiOD68uQsWhJ8rIiIiEisUcKPMzQ6KcLmcw5oLRfzpT06uXIlc7SIiIiLRQEsUosytHBQxnLkLFlTy9NNWLl6En/4UPv/5UfxQIiIiIqNIATcK3eygiOHOTUuD3FzYvJmw63hFRERExroxuURh7dq13HPPPUyePJmkpCRKSkqMLmlMeOgh+N//G1580ehKRERERCJnTAbcF154gQ8++ID29nZef/11nnrqKRobG40ua0z40pdAh5mJiIhILBuTSxSSk5P7fx0XF8fkyZOxWCwGViQiIiIi0WJMdnABCgoKSEhIYPHixbzyyitMnz7d6JJEREREJAqM2YBbXl5OV1cXbrcbh8PBe++9Z3RJMaG52YfDUUhqag5JSStITc3B4SjE5/MZXZqIiIjILYn6gFteXo7FYsFisZCTkzNgzGw2s3LlShYvXsyhQ4cMqjB2/NM/tXD//XmUlT2K13uYCxcO4fX+jLKyR8nMzLtpyPX5FI5FRETEeCMecDs7O3E6nWRnZ2O1WjGbzWzfvj3sazdv3ozNZiMhIYG5c+dSWVk54DUFBQV0dHTQ0dHBkSNHQr6P3+/nzjvvHOmPMu6cOrWH7u7dwEKuHfdrBhbS0LALp9MVdm5LSwuZmcMPxyIiIiIjZcQD7qVLlygpKaG7u5vc3FwATCZTyNeuWrWK/fv3s23bNo4fP86CBQvIz8/H7XaHff+LFy9SXV1NV1cXfr+fAwcO4PF4yMrKGumPMu5cuOAFMsKMZvCLX3jDzt2yZQ8NDcMLxyIiIiIjacR3UZg5cyaXL18GoLW1lX379oV83dGjRzl16hRut5u8vDwAlixZQlNTE4WFheTl5WE2h87fe/fuZcOGDZjNZtLS0jh06BA2m22kP8q44/fHcS2cXs/cNx5afb0XCBdiM6iv33mb1YmIiIjcmohuE9bb2xt27ODBg1gsFlavXj3gusPhYN26dXg8HjIzM2+Yd/fdd/PWW28NuZasrCwmTJiAzWbrD8P5+fnk5+cP+b1iVXx8D9BL6JAbICGhJ+zc2wnHIiIiMn643e7+n9Y3NzfT3NxMd3f3iN7DsH1wz549S3Jy8g1d2tmzZwNw7ty5kAF3uE6ePMm8efNG7P1iUXp6Cl6vh+Ayg+t5SE9PCTv3ZuE4OC4iIiLjXagG4+nTp5k/f/6I3cOwXRRaW1uZOnXqDdevXmttbR3tksY9l8uJ3f4sUAcE+q4GgDrs9iJcLmfYucHw6wkzOng4FhERERlJUb9NmIweq9VKXV0l69fXkJKynAcfXEFKynLWr6+hrq4S6yBn/A4Wjk2mIh57LHw4FhERERlJhi1RmDZtWsgubVtbW/+4jD6r1Upp6Z5hzaurq8TpdFFfvxO/P474+B7mzk0hEKgkNTV8OBYREREZSYYF3Dlz5uB2uwkEAgPW4Z45cwaAtLQ0o0qTYRpuOBYREREZSYYtUcjNzaWzs5Pq6uoB18vKyrDZbGRkhNuPVUREREQkvIh0cI8dO0ZXVxcdHR1AcEeEq0E2JyeHhIQEli5dSlZWFhs3bqS9vR273Y7b7ebEiROUl5eHPRxCRERERGQwEQm4mzZtoqmpCQieYlZVVUVVVRUmk4nGxkYSExMBqKmpoaioiK1bt9LW1kZycjIVFRWsWbMmEmWJiIiIyDgQkYDb2Nh4S6+bNGkSxcXFFBcXR6IMGQN6eiBOZ0CIiIjICNI2YWKYf/1XWLAALl0yuhIRERGJJQq4YpgvfxmamyE7Gy5fNroaERERiRUKuGKYL38ZTp2CpiZYuhR+9zsfDkchqak5JCWtIDU1B4ejEJ/PZ3SpIiIiMoYYtg+uCMDs2XDyJPzFX7SQmrqWTz/dDbgAExDA662ntjbvpiepiYiIiFylDq4Ybt48+OpX9/SF24UEwy0E//dcSEPDLpxO16Dv4fOp+ysiIiJB6uBKVPj9770EO7ehZFBfvzPs3JaWFhYtWktDg7q/IiIiog6uRAm/P45rndvrmfvGg3p64OWX4fBh+M1vYPPmPX3hdnjdXxEREYkt6uBKVIiP7wF6CR1yA33jQT4f/Jf/Egy6QcPv/oqIiEjsUQdXokJ6egrgCTPq6RsPuuce+NOf4IMPoK4OvvCFW+/+ioiISOxTwJWo4HI5sdufBeqAQN/VAFCH3V6Ey+Uc8Pq4OLDZYOFC+Pznr3Z/QxnY/RUREZHYpyUKEhWsVit1dZU4nS7q63fi98cRH99DenoKLtfgD4mlp6fg9XoIrsG93sDur4iIiMQ+BVyJGlarldLSPUOe53I5qa3No6FhF5BB8AcTAcDT1/2tHOFKRUREJJop4MqYdzvdXxEREYk9CrgSE4bb/RUREZHYo4fMZNzo6jK6AhERERkNCrgyLvh8kJICpaVGVyIiIiKRpoAr48L06bB0KXz72/DznxtdjYiIiESSAq6MCyYT/NM/wcMPw6OPgtdrdEUiIiISKQq4Mm7Ex8OBA5CYCH/5l/DHPxpdkYiIiESCAq6MK5Mnw5EjcOUKrFgBH39sdEUiIiIy0hRwZdz58z+Hw4fh3Dn4xjegRyf5ioiIxBQFXBmX5s2Digrw++FPfzK6GhERERlJYy7gnjt3joceeoi77roLu93Oq6++anRJMkYtXw5vvAETJ4LP58PhKCQ1NYekpBWkpubgcBTi8/mMLlNERESGaMydZPbEE0+wcuVK3nrrLX71q1+xZMkSFi9ezJe//GWjS5MxyGSClpYWFi1aS0PDbsAFmIAAXm89tbV51NXpuF8REZGxZMx1cM+fP09+fj4Ac+fOJTk5mXfffdfgqmQs27JlT1+4XUgw3ELwt8ZCGhp24XS6Bp2v7q+IiEh0GXMBNzs7m9deew2/34/H4+H9998nMzPT6LJkDKuv9wIZYUYz+OUvvfT2hh5taWkhMzOPsrJH8XoPc+HCIbzen1FW9iiZmXkKuSIiIgYYcwHX5XJRWlpKQkICX/nKV9ixYwczZswwuiwZw/z+OK51bq9n5sKFOO68E1JTg/vnbtoELhe8//7td39FRERk5EV1wC0vL8disWCxWMjJyeHjjz/m4Ycf5sUXX+TKlStcuHCB4uJiDh8+bHSpMobFx/cAYVq0BLjvvh5+8IPgKWgTJsAvfgG7dsH/+3837/4Gx0VERGQ0jWjA7ezsxOl0kp2djdVqxWw2s3379rCv3bx5MzabjYSEBObOnUtlZeWA1xQUFNDR0UFHRwdHjhzh3Llz+P1+Vq9ejclk4v7772f58uW8+eabI/kxZJxJT08BPGFGPTzySArPPAN79wZ3XfjNb+Cjj2DBgpt3f4PjIiIiMppGNOBeunSJkpISuru7yc3NBcBkCv2X/6pVq9i/fz/btm3j+PHjLFiwgPz8fNxud9j3t9vttLe3U1NTQ29vL01NTbzxxhvMmTNnJD+GjDMulxO7/VmgDgj0XQ0AddjtRbhczpDzTKabd3+D4yIiIjKaRnSbsJkzZ3L58mUAWltb2bdvX8jXHT16lFOnTuF2u8nLywNgyZIlNDU1UVhYSF5eHmbzjdl76tSpVFRU8Nxzz7F+/XosFgvr1q1jw4YNI/kxZJyxWq3U1VXidLqor9+J3x9HfHwP6ekpuFyDbxGWnp6C1+shuAb3ep6+7rCIiIiMpojtg9sb7rFz4ODBg1gsFlavXj3gusPhYN26dXg8nrA7Iyxbtoxly5aNaK0iVquV0tI9Q57ncjmprc2joWEXwbW4ZoLdXw+f+1wRu3dXDv4GIiIiMuIMOejh7NmzJCcn39ClnT17NhA8rWykt/7KyspiwoQJ2Gw2bDYbAPn5+f176ooMR7jub2JiCidOVPL3f29lz9Bzs4iISMxyu939S1Kbm5tpbm6mu7t7RO9hSMBtbW1l1qxZN1yfOnVq//hIO3nyJPPmzRvx9xUJ1/394Q9h82aYPx/WrjWgMBERkSgUqsF4+vRp5s+fP2L3iOptwkTGsqeegoIC+OY34f/8H6OrERERGT8MCbjTpk0L2aVta2vrHxcZ60wmeOUVSEqC3Fzo+99bREREIsyQgDtnzhzOnz9PIBAYcP3MmTMApKWlGVGWyIibOBEOHoQPP4QnnzS6GhERkfHBkICbm5tLZ2cn1dXVA66XlZVhs9nIyAh3MpTI2DNzZjDk7thhdCUiIiLjw4g/ZHbs2DG6urro6OgAgjsiXA2yOTk5JCQksHTpUrKysti4cSPt7e3Y7XbcbjcnTpygvLw87OEQImPVQw8ZXYGIiMj4MeIBd9OmTTQ1NQHBU8yqqqqoqqrCZDLR2NhIYmIiADU1NRQVFbF161ba2tpITk6moqKCNWvWjHRJIiIiIjKOjHjAbWxsvKXXTZo0ieLiYoqLi0e6BBEREREZx7RNmIiIiIjEFAVcEREREYkphpxkJiJB//Zv8Ic/+DhxwkV9vbf/qN/09BRcLidWq9XoEkVERMYcBVwRA/23/9bC4cNr6e3dDbgAExDA662ntjaPurpKhVwREZEh0hIFEQNZLHv6wu1CguEWgr8tF9LQsAun0zXofJ/Ph8NRSGpqDklJK0hNzcHhKMTn80W4chERkeilDq6IgX79ay/Bzm0oGdTX7ww7t6WlhUWL1tLQoO6viIjIZ6mDK2Igvz+Oa53b65n7xkPbsmVPX7gdXvdXREQkVingihgoPr4H6A0zGugbD+34cS8Q7ljrDOrrvbdZnYiIyNikgCtioPT0FMATZtTTNx7a7XR/RUREYpkCroiBXC4ndvuzQB0Q6LsaAOqw24twuZxh586YMfzur4iISCzTQ2YiBrJardTVVeJ0uqiv33ndPriDPySWnp6C1+shuAb3eoN3f0VERGKZAq6IwaxWK6Wle4Y8z+VyUlubR0PDLoJrcc0Eu78eEhOLcLkqR7hSERGRsUEBV2SMCtX9jYvr4Z57Uqis1BZhIiIyfingioxhw+3+ioiIxDI9ZCYiIiIiMUUBV0RERERiigKuiIiIiMQUBVwRERERiSkKuCLjwKVLcOSI0VWIiIiMDgVckXHgpZdg1Sr41a+MrkRERCTyFHBFxoHnn4e0NFizBtrbja5GREQkshRwRcaBz30OKivh4kX49reht9foikRERCJHAVdknJg1C/btCwbdV14xuhoREZHIicqA+/LLLzNv3jzuuOMOtm/fPmBs7dq13HPPPUyePJmkpCRKSkoMqlJk7FmzBjZuhL/6K/jNb4yuRkREJDKiMuDee++97Nixg5UrV2IymQaMvfDCC3zwwQe0t7fz+uuv89RTT9HY2GhQpSJjz9//PSQnw+rV0NFhdDUiIiIjLyoD7te//nWWLVvGlClT6L1usWBycjLx8fEAxMXFMXnyZCwWixFlioxJn/scHDgA8+bBH//ow+EoJDU1h6SkFaSm5uBwFOLz+YwuU0REZNjijS5gOAoKCqipqQGgoqKC6dOnG1yRyNjypS/B3r0tLFq0loaG3YALMAEBvN56amvzqKurxGq1GlypiIjI0EVlB/dmysvL6erqwu1243A4eO+994wuSWTM2bJlT1+4XUgw3ELwj4SFNDTswul0DTrf51P3V0REopPhAbe8vByLxYLFYiEnJ+eW55nNZlauXMnixYs5dOhQBCsUiU319V4gI8xoRt94aC0tLWRm5lFW9ihe72EuXDiE1/szysoeJTMzTyFXREQMNeSA29nZidPpJDs7G6vVitlsvmGng8++dvPmzdhsNhISEpg7dy6VlZUDXlNQUEBHRwcdHR0cCXGW6PUPmV3P7/dz5513DvVjiIx7fn8c1zq31zNz4UIcWVnwne+AywU//WnwJLT29tvv/oqIiETSkAPupUuXKCkpobu7m9zcXCB8CF21ahX79+9n27ZtHD9+nAULFpCfn4/b7R70Hj09PXz66af4/X66u7v59NNPCQQCXLx4kerqarq6uvD7/Rw4cACPx0NWVtZQP4bIuBcf3wOEO/EhwJQpPUyZAh4P7NwJ//k/Bx9MmzIFDh8efvdXREQk0ob8kNnMmTO5fPkyAK2trezbty/k644ePcqpU6dwu93k5eUBsGTJEpqamigsLCQvLw+zOXS+/sEPfsCOHTv6/3vXrl2UlZXxta99jb1797JhwwbMZjNpaWkcOnQIm8021I8hMu6lp6fg9XoIdmGv52H58hRKS4P/1dsLra3wu98F/9myZfDub7A7LCIiYozb2kXh+i28PuvgwYNYLBZWr1494LrD4WDdunV4PB4yMzNDzt22bRvbtm0LOfbWW28Nq9asrCwmTJiAzWbrD8T5+fnk5+cP6/1ExjqXy0ltbR4NDbsIdmPNQADwYLcX4XJdW05kMsH06cF/0tPhBz+42v0NFXIDfd1hERGRG7nd7v6f5jc3N9Pc3Ex3d/eI3iNi24SdPXuW5OTkG7q0s2fPBuDcuXNhA24knDx5knnz5o3a/USindVqpa6uEqfTRX39Tvz+OOLje0hPT8HlGnyLsJt1f9PTUyJWt4iIjG2hGoynT59m/vz5I3aPiAXc1tZWZs2adcP1qVOn9o+LiLGsViulpXuGPG+w7u/MmQO7vyIiIqNtTB70ICLGCtf9veOOFNrbK5k4UQdEiIiIcSIWcKdNmxayS9vW1tY/LiJjV6ju77vvBndaeOYZ+PGPDSpMRETGvYgd9DBnzhzOnz9PIBAYcP3MmTMApKWlRerWImKQpCQoLoZXXoF//VejqxERkfEqYgE3NzeXzs5OqqurB1wvKyvDZrORkRFuD00RGcs2bICVK4P//sMfjK5GRETGo2EtUTh27BhdXV10dHQAwR0RrgbZnJwcEhISWLp0KVlZWWzcuJH29nbsdjtut5sTJ05QXl5+0xPKRGRsMpmgpATmzIH16+H4cQiz5bWIiEhEDCvgbtq0iaamJiB4illVVRVVVVWYTCYaGxtJTEwEoKamhqKiIrZu3UpbWxvJyclUVFSwZs2akfsEIhJ1pk+Hn/wEsrPhhz+Ep582uiIRERlPhhVwGxsbb+l1kyZNori4mOLi4uHcRkTGsKys4MNmL74ITz4JEycaXZGIiIwX+sGhiETM7t3w7/+ucCsiIqNLAVdEIubP/gzuvdfoKkREZLxRwBURERGRmKKAKyIiIiIxRQFXRERERGKKAq6IiIiIxBQFXBEZdW+8AefOGV2FiIjEqmHtgysiMlx+PxQVQW+vj/nzXbz9the/P474+B7S01NwuZxYrVajyxQRkTFMAVdERlV8PPzDP7TwH//jWrze3YALMAEBvN56amvzqKurVMgVEZFh0xIFERl1+/fvAXYDCwmGWwj+cbSQhoZdOJ2uQef7fD4cjkJSU3NISlpBamoODkchPp8vsoWLiMiYoA6uiIy6+novwc5tKBnU1+8MO7elpYVFi9bS0KDur4iIhKYOroiMOr8/jmud2+uZ+8Zv1NMDhYV7+sLt8Lq/IiIS+xRwRWTUxcf3AL1hRgN94zd6+23Yv98LZISZm9HXHRYRkfFMAVdERl16egrgCTPq6Ru/0Z//OcyYMbzur4iIjB8KuCIy6lwuJ3b7s0AdEOi7GgDqsNuLcLmcIed94Qswffrwur8iIjJ+6CEzERl1VquVurpKnE4X9fU7r9sHd/CHxNLTU/B6PQTX4F4vfPdXRETGDwVcETGE1WqltHTPkOe5XE5qa/NoaNhFcC2umWD319PX/a0c4UpFRGSsUcAVkTHldrq/IiIyPijgisiYM9zur4iIjA96yExEREREYooCroiIiIjEFAVcEREREYkpURlwX375ZebNm8cdd9zB9u3bB4ydO3eOhx56iLvuugu73c6rr75qUJUiIiIiEo2iMuDee++97Nixg5UrV2IyDTyx6IknniA7O5sPP/yQ6upqnn76ad555x2DKhWRaFVfD1/7Gnz8sdGViIjIaIvKgPv1r3+dZcuWMWXKFHp7B55YdP78efLz8wGYO3cuycnJvPvuu0aUKSJR7K67oLYW/uZvjK5ERERGW1QG3MFkZ2fz2muv4ff78Xg8vP/++2RmZhpdlohEmQcfhL/9W/iHf4Cf/9zoakREZDSNuYDrcrkoLS0lISGBr3zlK+zYsYMZM2YYXZaIRKHvfhf+w3+Ab34TPvrI6GpERGS0GB5wy8vLsVgsWCwWcnJyBn1tV1cXDz/8MC+++CJXrlzhwoULFBcXc/jw4VGqVkTGErMZSkvh8mV4+mmjqxERkdEy5IDb2dmJ0+kkOzsbq9WK2Wy+YaeDz7528+bN2Gw2EhISmDt3LpWVA8+JLygooKOjg46ODo4cOXLDe3z2ITOv14vf72f16tWYTCbuv/9+li9fzptvvjnUjyEi48QXvwjFxcGg+7OfGV2NiIiMhiEH3EuXLlFSUkJ3dze5ubkAN+x0cNWqVavYv38/27Zt4/jx4yxYsID8/Hzcbveg9+jp6eHTTz/F7/fT3d3Np59+SiAQwG63097eTk1NDb29vTQ1NfHGG28wZ86coX4MERlHHA5Ytgy+9S1oazO6GhERibT4oU6YOXMmly9fBqC1tZV9+/aFfN3Ro0c5deoUbrebvLw8AJYsWUJTUxOFhYXk5eVhNofO1z/4wQ/YsWNH/3/v2rWLsrIyHn/8cSoqKnjuuedYv349FouFdevWsWHDhqF+DBEZR0wmKCmBqqrg7goiIhLbTL3X78M1BJcuXWLGjBls27aNrVu3Dhj71re+xYEDB7h8+fKAIFtRUcG6dev4xS9+MSq7H5w+fZr58+czdepUJkyYgM1mw2azAZCfn9+/5ZiIiIiIRJ7b7e7/aX5zczPNzc10d3fT1tbG22+/zbx58277HkPu4N6qs2fPkpycfEOXdvbs2UDwRLLR3N7r5MmTI/KFicjY5/P5cDpd1Nd78fvjiI/vIT09BZfLidVqNbo8EZGYFqrBeLUhOVIiFnBbW1uZNWvWDdenTp3aPy4iMtpaWlpYtGgtDQ27ARdgAgJ4vfXU1uZRV1epkCsiMsYZvk2YiMho2rJlT1+4XUgw3ELwj8KFNDTswul0DTrf5/PhcBSSmppDUtIKUlNzcDgK8fl8Ea5cRERuVcQ6uNOmTQvZpW3re4R52rRpkbq1iEhY9fVegp3bUDL45S930tsbfDDteur+ioiMDRHr4M6ZM4fz588TCAQGXD9z5gwAaWlpkbq1iEhYfn8c1zq31zNz4UIcU6bAvHmwZg08+yycPx8cvd3ur4iIjI6IBdzc3Fw6Ozuprq4ecL2srAybzUZGRkakbi0iElZ8fA8QbvOYAPfd18Pzz8OCBdDaCuXl0NwcHA12f8P92ZXRNy4iIkYb1hKFY8eO0dXVRUdHBxDcEeFqkM3JySEhIYGlS5eSlZXFxo0baW9vx26343a7OXHiBOXl5WEPhxARiaT09BS8Xg/BLuz1PDzySAqFhaHn3qz7GxwXERGjDSvgbtq0iaamJiB4illVVRVVVVWYTCYaGxtJTEwEoKamhqKiIrZu3UpbWxvJyclUVFSwZs2akfsEIiJD4HI5qa3No6FhF8FurBkIAB7s9iJcrsqwc691f0OF3EDfuIiIGG1YAbexsfGWXjdp0iSKi4spLi4ezm1EREac1Wqlrq6ybx/cndftgzv4Q2I36/6mp6dErG4REbl1EdtFQUQkWlmtVkpL9wx53mDd38TEwbu/IiIyehRwRURuUajur9ncw/vvp3DffZVMm6YtwkREooECrojIEITq/v785/DII/DjH8PGjQYVJiIi/XSSmYjIbXr4Yfj2t8HphL7nb0VExEAKuCIiI8Dlgrvugm99C3rDbbMrIiKjQgFXRGQETJkCr7wCJ09CRYXR1YiIjG9agysiMkL+03+CykpYscLoSkRExjcFXBGREaRzbEREjKclCiIiIiISUxRwRUR5dUghAAAgAElEQVRERCSmKOCKiIiISExRwBURERGRmKKAKyISYT090NlpdBUiIuOHAq6ISIStXBk86UxEREaHAq6ISITl54PbDW+8YXQlIiLjgwKuiEiE5efD8uXwne9AW5vR1YiIxD4FXBGRCDOZ4Ec/gk8+gWeeMboaEZHYp5PMRERGwb33wksvwTe/6aO52cUf/uDF748jPr6H9PQUXC4nVqvV6DJFRGKCAq6IyCj5y79sISFhLadO7QZcgAkI4PXWU1ubR11dpUKuiMgI0BIFEZFR8jd/s4dPPtkNLCQYbiH4x/BCGhp24XS6Bp3v8/lwOApJTc0hKWkFqak5OByF+Hy+CFcuIjK2qIMrIjJK6uu9BDu3oWRQX78z7NyWlhYWLVpLQ4O6vyIiN6OAKyIySvz+OK51bq9n7hsP+ou/CD6UNmsWfOlL8D/+x56+cLtwwJzPdn9LS/dErHYRkbEk6pYoXLlyBYfDQWJiIlOmTCEzM5O6urr+8bVr13LPPfcwefJkkpKSKCkpMbBaEZFbFx/fA/SGGQ30jQetWAEpKfDee8EdGGprvUBGmLkZfd1hERGBKAy4fr+fBx54gF/+8pd89NFHbNy4kRUrVvDpp58C8MILL/DBBx/Q3t7O66+/zlNPPUVjY6PBVYuI3Fx6egrgCTPq6RsPeuYZKC2F2lr44x9h1qxb7/6KiIx3URdwJ06cyPPPP899990HwOOPP04gEOC3v/0tAMnJycTHB1dWxMXFMXnyZCwWi2H1iojcKpfLid3+LFAHBPquBoA67PYiXC5n2Ll33HHr3V8RkfEu6tfgvvPOO3zyySfY7fb+awUFBdTU1ABQUVHB9OnTjSpPROSWWa1W6uoqcTpd1NfvvG4f3MEfEktPT8Hr9TBwDe5VA7u/IiLjXVQH3I8//pjHHnuM559/nokTJ/ZfLy8vJxAIcOjQIRwOB7/+9a9JTEw0sFIRkVtjtVqH9TCYy+WktjaPhoZdBNfimgl2fz193d/KEa5URGTsMnyJQnl5ORaLBYvFQk5OTv/17u5uVq9eTVpaGt///vdvmGc2m1m5ciWLFy/m0KFDo1myiMiou9r9Xb++hpSU5Tz44ApSUpazfn2NtggTEbnOkANuZ2cnTqeT7OxsrFYrZrOZ7du3h33t5s2bsdlsJCQkMHfuXCorB3YZCgoK6OjooKOjgyNHjgAQCAR47LHHuOOOO3j11VcHrcfv93PnnXcO9WOIiIw5V7u/584d4d13D3Hu3BFKS/co3IqIXGfIAffSpUuUlJTQ3d1Nbm4uACZT6Cd7V61axf79+9m2bRvHjx9nwYIF5Ofn43a7B73Hk08+ycWLF6moqMBsvlbixYsXqa6upqurC7/fz4EDB/B4PGRlZQ31Y4iIxJxPPjG6AhGR6DDkNbgzZ87k8uXLALS2trJv376Qrzt69CinTp3C7XaTl5cHwJIlS2hqaqKwsJC8vLwB4fWqpqYmXn31VRISEgY8PHb8+HFmzZrF3r172bBhA2azmbS0NA4dOoTNZhvqxxARiSmbN4PXC2++CWF6DiIi48ZtPWTW2xtuyxo4ePAgFouF1atXD7jucDhYt24dHo+HzMzMG+Z98YtfJBAI3HD9qrfeemv4BYuIxKilS+GHP4TXXoPHHze6GhERY0VsF4WzZ8+SnJx8Q5d29uzZAJw7dy5kwI2UrKwsJkyYgM1m6+/45ufnk5+fP2o1iIhEytKlUFAATz8d/PWMGUZXJCISmtvt7l+u2tzcTHNzM93d3SN6j4gF3NbWVmbNmnXD9alTp/aPj6aTJ08yb968Ub2niMhoeuklOH48GHLLy42uRkQktFANxtOnTzN//vwRu4fh24SJiMjIsFqDIfdf/gWOHTO6GhER40Qs4E6bNi1kl7atra1/XERERtY3vgFZWfCd70Bnp9HViIgYI2IBd86cOZw/f/6GB8bOnDkDQFpaWqRuLSIybplM8KMfgc8HO3YYXY2IiDEiFnBzc3Pp7Oykurp6wPWysjJsNhsZGRmRurWIyLj2wANw4AA884zRlYiIGGNYD5kdO3aMrq4uOjo6gOCOCFeDbE5ODgkJCSxdupSsrCw2btxIe3s7drsdt9vNiRMnKC8vD3s4hIiI3L5ly4yuQETEOMMKuJs2baKpqQkInmJWVVVFVVUVJpOJxsZGEhMTAaipqaGoqIitW7fS1tZGcnIyFRUVrFmzZuQ+gYiIiIjIZwwr4DY2Nt7S6yZNmkRxcTHFxcXDuY2IiIiIyJBpmzARERERiSkKuCIiIiISUyJ2kpmIiEQXn8/HX/+1i7ff9uL3xxEf30N6egoulxOr1Wp0eSIiI0YBV0RkHGhpaeHBB9fy0Ue7ARdgAgJ4vfXU1uZRV1epkCsiMUNLFERExoEtW/b0hduFBMMtBP8KWEhDwy6cTteg830+Hw5HIampOSQlrSA1NQeHoxCfzxfhykVEhk4dXBGRcaC+3kuwcxtKBvX1O8PObWlpYdGitTQ0qPsrImODOrgiIuOA3x/Htc7t9cy0tsaFnbtly56+cDu87q+IyGhTwBURGQfi43uA3jCjAe68syfs3GD3N9zx6hl94yIi0UMBV0RkHEhPTwE8YUY9fPWrKWHn3qz7GxwXEYkeCrgiIuOAy+XEbn8WqAMCfVcDQB12exEulzPs3Jt1f4PjIiLRQw+ZiYiMA1arlbq6SpxOF/X1O6/bB3fwh8TS01Pwej0E1+Bez9PXHRYRiR4KuCIi44TVaqW0dM+Q57lcTmpr82ho2EVwLa6ZYPfXwwMPFOFyVY5wpSIit0cBV0REBhWq+3vlSg+//30KTz6pLcJEJPoo4IqIyE2F6v7u2QOPPGJQQSIig1DAFRGRYSksNLoCEZHQtIuCiIiIiMQUBVwRERERiSkKuCIiIiISUxRwRURERCSmKOCKiIiISExRwBURkRHR2wsvvAA/+YnRlYjIeKeAKyIiI8JkgsZGeOYZ8PmMrkZExjMFXBERGTEvvhjs5H7ve0ZXIiLjWdQF3CtXruBwOEhMTGTKlClkZmZSV1fXP37u3Dkeeugh7rrrLux2O6+++qqB1YqIyGfNmAEuF+zfD//zfxpdjYiMV1EXcP1+Pw888AC//OUv+eijj9i4cSMrVqzg008/BeCJJ54gOzubDz/8kOrqap5++mneeecdg6sWEZGrvvlN+MpX4DvfgT/9yehqRGQ8irqAO3HiRJ5//nnuu+8+AB5//HECgQC//e1vATh//jz5+fkAzJ07l+TkZN59913D6hURkYHMZvjRj+B3v4O/+zujqxGR8SjqAu713nnnHT755BPsdjsA2dnZvPbaa/j9fjweD++//z6ZmZkGVykiIp+VmgqFhbBrF1y4YHQ1IjLexBtdwGA+/vhjHnvsMZ5//nkmTpwIgMvl4pFHHmHXrl0AvPzyy8yYMcPIMkVEJITnnoPz56FvhZmIyKgxvINbXl6OxWLBYrGQk5PTf727u5vVq1eTlpbG97//fQC6urp4+OGHefHFF7ly5QoXLlyguLiYw4cPG1W+iIiEMXEiHDwIc+YYXYmIjDdDDridnZ04nU6ys7OxWq2YzWa2b98e9rWbN2/GZrORkJDA3LlzqaysHPCagoICOjo66Ojo4MiRIwAEAgEee+wx7rjjjgG7JHi9Xvx+P6tXr8ZkMnH//fezfPly3nzzzaF+DBERGUU+nw+Ho5DU1BySklaQmpqDw1GITxvmikgEDDngXrp0iZKSErq7u8nNzQXAZDKFfO2qVavYv38/27Zt4/jx4yxYsID8/Hzcbveg93jyySe5ePEiFRUVmM3XSrTb7bS3t1NTU0Nvby9NTU288cYbzFF7QEQkarW0tJCZmUdZ2aN4vYe5cOEQXu/PKCt7lMzMPIVcERlxQ16DO3PmTC5fvgxAa2sr+/btC/m6o0ePcurUKdxuN3l5eQAsWbKEpqYmCgsLycvLGxBer2pqauLVV18lISGB6dOn918/fvw4ixcvpqKigueee47169djsVhYt24dGzZsGOrHEBGRUbJlyx4aGnYDCz9z1QwspKFhF06ni9LSPQZVJyKx6LYeMuvt7Q07dvDgQSwWC6tXrx5w3eFwsG7dOjweT8jdD774xS8SCATCvu+yZctYtmzZkGvNyspiwoQJ2Gw2bDYbAPn5+f1bjomISGTU13sBV5jRDOrrdw463+fz4XS6qK/34vfHER/fQ3p6Ci6XE6vVOuL1ikhkud3u/p/mNzc309zcTHd394jeI2K7KJw9e5bk5OQburSzZ88GgieSjeb2XidPnmTevHmjdj8REQny++OA0EvZwExTUxxOJ3zpS/Dgg8F/f+ELYDIFlzcsWrS2rwPs6nufAF5vPbW1edTVVSrkiowxoRqMp0+fZv78+SN2j4jtotDa2srUqVNvuH71Wmtra6RuLSIiUSQ+vgcI9xO/ANBDdXXw5LO/+Auw2aCiIjg6cHnD1ZA8cHmDiMj1DN8mTEREYlt6egrgCTPqYfXqFH73O/j4Y/B64Y03gkEXri5vyAgzN6NvXERkoIgF3GnTpoXs0ra1tfWPi4hI7HO5nNjtzwJ1BDu29P27Dru9CJfLCcCf/RkkJ8OKFcElCnDz5Q3BcRGRgSIWcOfMmcP58+dveGDszJkzAKSlpUXq1iIiEkWsVit1dZWsX19DSspyHnxwBSkpy1m/vuama2hvtrwhOC4iMlDEHjLLzc2lpKSE6upq1qxZ03+9rKwMm81GRka4HzmJiEissVqtw9oKLD09Ba/Xw8Atxq7y9C1/EBEZaFgB99ixY3R1ddHR0QEEd0Sorq4GICcnh4SEBJYuXUpWVhYbN26kvb0du92O2+3mxIkTlJeXhz0cQkRE5CqXy0ltbR4NDbsIrsU1E1ze4MFmK8Llqhz8DURkXBpWwN20aRNNTU1A8BSzqqoqqqqqMJlMNDY2kpiYCEBNTQ1FRUVs3bqVtrY2kpOTqaioGNDRFRERCefq8obgPrg78fvjiIvroaUlhSlTKvn857VFmIjcaFgBt7Gx8ZZeN2nSJIqLiykuLh7ObUREREIub3j7bViwAP7xH2HzZoMKE5GopW3CRERkzJk/HzZuhK1b4Q9/MLoaEYk2CrgiIjIm7dwJCQnw139tdCUiEm0UcEVEZEz6/Odhzx74+c/B5zO6GhGJJgq4IiIyZj32GFy4AINspSsi45ACroiIjFkmE9x1l9FViEi0UcAVERERkZiigCsiIiIiMUUBV0RERERiigKuiIiIiMQUBVwREYk5Fy4YXYGIGEkBV0REYsprr8Hs2fDb3xpdiYgYRQFXRERiyqOPwr33wne/C729RlcjIkZQwBURkZgycSLs3QsnTsBPf2p0NSJiBAVcERGJOcuXw4oVsHkzdHYaXY2IjDYFXBERiUk//CG0tcGOHUZXIiKjLd7oAkRERCJh5kwoKoIXXvBx4YKL3/7Wi98fR3x8D+npKbhcTqxWq9FlikgEqIMrIiIx6/HHWzCb83jjjUfxeg9z4cIhvN6fUVb2KJmZefh8PqNLFJEIUMAVEZGYtXXrHrq7dwMLAVPfVTOwkIaGXTidrkHn+3w+HI5CUlNzSEpaQWpqDg5HoYKxSJTTEgUREYlZ9fVeIFyIzaC+fmfYuS0tLSxatJaGht1972ECAni99dTW5lFXV6klDiJRSh1cERGJWX5/HNc6t9cz4/fH8dOfBvfM3bsXjh2Dhgbw+2HLlj194XZ43V8RMY46uCIiErPi43uAXkKH3ADx8T34fPBv/wb79sGf/hQcmTABTKbhd39FxFjq4IqISMxKT08BPGFGPaSnp/Cd78CZM/Dxx9DUBCdPQnExTJx48+6viESnqAy4a9eu5Z577mHy5MkkJSVRUlJyS2MiIiKf5XI5sdufBeqAQN/VAFCH3V6Ey+Xsf63ZDImJ8MgjsGkT3Hvv1e5vKIG+7rCIRKOoDLgvvPACH3zwAe3t7bz++us89dRTNDY23nRMRETks6xWK3V1laxfX0NKynIefHAFKSnLWb++5qYPid1K91dEolNUrsFNTk7u/3VcXByTJ0/GYrHcdExEROR6VquV0tI9Q57ncjmprc2joWEXkEGwJxQAPH3d38oRrlRERkpUBlyAgoICampqAKioqGD69Om3NCYiIjISrnZ/nU4X9fU7rzsFTVuEiUSzqA245eXlBAIBDh06hMPh4Ne//jWJiYk3HRMRERkpw+3+ioixDF+DW15ejsViwWKxkJOTM2DMbDazcuVKFi9ezKFDh255TEREJJK2boWXXjK6ChEJZ8gBt7OzE6fTSXZ2NlarFbPZzPbt28O+dvPmzdhsNhISEpg7dy6VlQPXLBUUFNDR0UFHRwdHjhwJ+T5+v58777xzyGMiIiKR0N4Ozz8Pzc1GVyIioQw54F66dImSkhK6u7vJzc0FwGQKvU/gqlWr2L9/P9u2beP48eMsWLCA/Px83G532Pe/ePEi1dXVdHV14ff7OXDgAB6Ph6ysrEHHRERERsu2bTBpEhQWGl2JiIQy5DW4M2fO5PLlywC0trayb9++kK87evQop06dwu12k5eXB8CSJUtoamqisLCQvLw8zObQ+Xrv3r1s2LABs9lMWloahw4dwmazcfHixbBjIiIio+Wuu+Dv/g4cDnjySViyxOiKROSzbushs97ecBtgw8GDB7FYLKxevXrAdYfDwbp16/B4PGRmZt4w7+677+att94K+Z6DjYmIiIymxx+HH/8Yvvtd+NWvID5qH9sWGX8i9tvx7NmzJCcn39ClnT17NgDnzp0LGXAjJSsriwkTJmCz2fo7vvn5+eTn549aDSIiEjvMZvjHf4QFC+C//3d46imjKxIZG9xud/9y1ebmZpqbm+nu7h7Re0Qs4La2tjJr1qwbrk+dOrV/fDSdPHmSefPmjeo9RUQkts2fD9/+dvCBs7w8uPtuoysSiX6hGoynT59m/vz5I3YP/UBFRETkNuzaBb//PXz4oQKuSLSIWMCdNm1ayC5tW1tb/7iIiMhYN20aHD9udBUi8lkRO+hhzpw5nD9/nkAgMOD6mTNnAEhLS4vUrUVERERkHItYwM3NzaWzs5Pq6uoB18vKyrDZbGRkZETq1iIiIiIyjg1ricKxY8fo6uqio6MDCO6IcDXI5uTkkJCQwNKlS8nKymLjxo20t7djt9txu92cOHGC8vLysIdDiIiIiIjcjmEF3E2bNtHU1AQETzGrqqqiqqoKk8lEY2MjiYmJANTU1FBUVMTWrVtpa2sjOTmZiooK1qxZM3KfQERERETkM4YVcBsbG2/pdZMmTaK4uJji4uLh3EZEREREZMgitgZXRERkvPr0U9i3D657zlpERokCroiIyAj793+Hb30LXn/d6EpExicd9CAiIjLCvvpVWLsWvvc9H2++6eLXv/bi98cRH99DenoKLpcTq9VqdJkiMUsBV0REJAK2bGmhsnIt//IvuwEXYAICeL311NbmUVdXqZArEiFaoiAiIhIBP/zhHnp7dwMLCYZbCP61u5CGhl04nS7jihOJcQq4IiIiEVBf7wXCHWqU0Tcens/nw+EoJDU1h6SkFaSm5uBwFOLz+Ua8VpFYoyUKIiIiEeD3x3Gtc3s9c994aC0tLSxatJaGBi1vEBkOdXBFREQiID6+B+gNMxroGx9o3z548knIydnTF261vEFkOBRwRUREIiA9PQXwhBn19I0P1N4O9fXw9tu3t7xBZLxTwBUREYkAl8uJ3f4sUAdcPfEhANRhtxfhcjlvmPPMM/CrX8GXvjT85Q0iojW4IiIiEWG1Wqmrq8TpdFFfv/O6fXAHX0N7bXlDqJAbenmDiFyjgCsiIhIhVquV0tI9Q56Xnp6C1+shuAb3eqGXN4jINVqiICIiEmWGs7xBRK5RB1dERCTK3M7yBhFRwBUREYlKw13eICJaoiAiIiIiMUYBV0RERERiigKuiIjIGBUI3Pw1IuORAq6IiMgY09sLq1fD7t1GVyISnRRwRURExhiTCWbOhP/6X+GDD4yuRiT6KOCKiIiMQc8/DxYLOLUlrsgNFHBFRETGoMmTgx1ctxv+1/8yuhqR6KKAKyIiMkY98QQsWABPPQU9PUZXIxI9ojLgrl27lnvuuYfJkyeTlJRESUlJ/9i5c+d46KGHuOuuu7Db7bz66qsGVioiImIcsxn27oVf/Qr++Z+NrkYkekRlwH3hhRf44IMPaG9v5/XXX+epp56isbERgCeeeILs7Gw+/PBDqqurefrpp3nnnXcMrlhERMQYCxfCY4/Bs8/Chx8aXY1IdIjKo3qTk5P7fx0XF8fkyZOxWCwAnD9/nsrKSgDmzp1LcnIy7777Ll/+8pcNqVVERMRof/u38LnPgd9vdCUi0SEqO7gABQUFJCQksHjxYl555RWmT58OQHZ2Nq+99hp+vx+Px8P7779PZmamwdWKiIgY59574ZVXoO+vSpFxL2oDbnl5OV1dXbjdbhwOB++99x4ALpeL0tJSEhIS+MpXvsKOHTuYMWOGwdWKiIiISLQwPOCWl5djsViwWCzk5OQMGDObzaxcuZLFixdz6NAhurq6ePjhh3nxxRe5cuUKFy5coLi4mMOHDxtUvYiIiIhEmyEH3M7OTpxOJ9nZ2VitVsxmM9u3bw/72s2bN2Oz2UhISGDu3Ln962evKigooKOjg46ODo4cORLyffx+P3feeSderxe/38/q1asxmUzcf//9LF++nDfffHOoH0NEREREYtSQA+6lS5coKSmhu7ub3NxcAEwmU8jXrlq1iv3797Nt2zaOHz/OggULyM/Px+12h33/ixcvUl1dTVdXF36/nwMHDuDxeMjKysJut9Pe3k5NTQ29vb00NTXxxhtvMGfOnKF+DBERERGJUUPeRWHmzJlcvnwZgNbWVvbt2xfydUePHuXUqVO43W7y8vIAWLJkCU1NTRQWFpKXl4fZHDpf7927lw0bNmA2m0lLS+PQoUPYbDYAKioqeO6551i/fj0Wi4V169axYcOGoX4MERGRmObz+XA6XdTXe/H744iP7yE9PQWXy4nVajW6PJGIuq1twnp7e8OOHTx4EIvFwurVqwdcdzgcrFu3Do/HE3L3g7vvvpu33nor7PsuW7aMZcuWDbnWrKwsJkyYgM1m6w/L+fn55OfnD/m9REREollLSwvp6WtpatoNuAATEMDrrae2No+6ukqFXDGM2+3u/2l+c3Mzzc3NdHd3j+g9IrYP7tmzZ0lOTr6hSzt79mwgeCLZaG7vdfLkSebNmzdq9xMRETHKli17+sLtws9cNQMLaWjYhdPporR0j0HVyXgXqsF4+vRp5s+fP2L3iNguCq2trUydOvWG61evtba2RurWIiIi41p9vRfICDOa0Tcens/nw+EoJDU1h6SkFaSm5uBwFOLz+Ua8VpFIiMqTzERERGT4/P44gssSQjFz6VJc2LktLS0sWrSWhgYtb5CxK2Id3GnTpoXs0ra1tfWPi4iIyMiLj+8Bwj0nE6C3tyfs3C1b9vSF24VcC8kDlzeIRLuIBdw5c+Zw/vx5AoHAgOtnzpwBIC0tLVK3FhERGdfS01MAT5hRDzk5KWHn3u7yBpFoELGAm5ubS2dnJ9XV1QOul5WVYbPZyMgI95tHREREbofL5cRufxaoA642mgJAHXZ7ES6XM+zcmy1vCI6LRLdhrcE9duwYXV1ddHR0AMEdEa4G2ZycHBISEli6dClZWVls3LiR9vZ27HY7brebEydOUF5eHvZwCBEREbk9VquVurrKvn1wd163D+7ga2ivLW8I9fd0oG9cJLoNK+Bu2rSJpqYmIHiKWVVVFVVVVZhMJhobG0lMTASgpqaGoqIitm7dSltbG8nJyVRUVLBmzZqR+wQiIiJyA6vVOqytwNLTU/B6PQzcYuwqT9/yB5HoNqyA29jYeEuvmzRpEsXFxRQXFw/nNiIiIjLKXC4ntbV5NDTsIrgW10xweYOHmTOLcLkqjS1Q5BZomzARERHpF2p5A/Twu9+lsG6dtgiTsUEBV0RERAYItbzh9GmYM8eggkSGSAFXREREbkqn3ctYErFtwkREREREjKCAKyIiIiIxRQFXRERERGKKAq6IiIiIxBQFXBERERmW3l6jKxAJTQFXREREhuz//l+YOxfeecfoSkRupIArIiIiQ3bfffDRR/C97/3/9u4/quoq3//48xwwB+n4+2BFmYZp4I8ltkTJMpuUnEt6BTNCbiu4uVqjNSWVaJo/0/QebTTX7a4aaXCZzIFgcI35I9PGGpu4MGV6VUyNITK7VxBQfnwtQc73jwMoAgp44HM4vB5rnWXs/fmcz/vstW29fbPP3kZHItKQElwRERFpsV/9CtauhZ07Yc8eo6MRqU8JroiIiLTK9OkwfjzEx0NVldHRiFyhBFdERERaxWSC9eud63DffdfoaESuUIIrIiIirTZqFMTFwdKlUFxsdDQiTkpwRURE5KasWgWXLsGKFUZHIuLkbXQAIiIi0rHddpszya2sNDoSEScluCIiInLTXnzR6AhErtASBRERERHxKEpwRURERMSjaImCiIiIuFRhYSEJCTays3OoqvLC2/syISFB2GwJWK1Wo8OTTkAJroiIiLhMQUEBDzzwFLm5bwI2wARUk5OTzYEDUWRmpirJlTanJQoiIiLiMvPnr61JbsfiTG7BmW6MJTd3FQkJtuveX1hYSFzcPIYODWfIkKkMHRpOXNw8CgsL2zhy8SRuneBmZmZiNptZtWpVXdtTTz3FbbfdRvfu3RkyZAibNm0yMEIRERG5WnZ2DjCmid4xfPFFTpPH+hYUFBAaGsXmzdPJydnByZPbycn5iM2bpxMaGqUkV5rNbRPc6upq4uPjCQ0NxWQy1bUvXbqUH3/8kdLSUrZu3cqLL75IXl6egZGKiIhIraoqL65Ubq9l5rvvvPDxgfvug1deqd97s9VfkVpum+C+996SeZ4AACAASURBVN57jBs3jsGDB+NwOOraAwMD8fZ2Lh328vKie/fuWCwWo8IUERGRq3h7XwYcTfRWc/fdl3n7bZg8Gfr1q997o+qvs1/kxtzyS2ZFRUVs3LiRrKwsXnrppQb9MTExZGRkAJCSkkLfvn3bO0QRERFpREhIEDk5WTirsNfK4pFHgpgzp/F7b1T9dfaL3JhbVnBfe+01XnnlFbp37w5Qb4kCQHJyMhUVFdjtduLi4vjhhx+MCFNERESuYbMlEBCwEMgEqmtaq4FMAgIWYbMlNHnvjaq/zn6RGzM8wU1OTsZisWCxWAgPD+frr7/mm2++4dlnnwXA4XDUW6JQy2w2M23aNMaNG8f27dvbO2wRERFphNVqJTMzldjYDIKCpjB48FSCgqYQG5txwy3CQkKCgKwmerNq+kVurMVLFMrLy1mxYgWHDh3im2++oaioiKVLl7J06dJGr3399ddJS0ujuLiY++67jwULFhAVFVV3TUxMDDExMXU/v/322+Tk5ODn51f3Hl5eXpw6dYrNmzc3eEZVVRW33nprSz+GiIiItBGr1UpS0toW32ezJXDgQBS5uatwrsU146z+ZtVUf1NdHKl4qhZXcM+dO8emTZuorKwkIiICaLiEoFZkZCRbtmxh2bJlfPzxx4wePZro6GjsdnuT7z9r1ixOnTrF4cOHOXToEFOnTuWFF15g/fr1nD17lvT0dCoqKqiqquLDDz8kKyuLSZMmtfRjiIiIiJtprPp7551TgAxWr9YBEdJ8La7gDhgwgJKSEsD5ZbDExMRGr9u1axf79u3DbrfXVWwffvhh8vPzmTdvHlFRUZjNDfNrX19ffH19637u1q0b3bt3p1evXpw9e5aNGzcya9YszGYzw4YNY/v27fj7+7f0Y4iIiIgburb663DAQw/BG29AZCR46Xtm0gw3tYtCY2tja23btg2LxcKMGTPqtcfFxTFz5kyysrIIDQ294TOSkpLq/rtfv3787W9/a1WskyZNokuXLvj7+9clxNHR0URHR7fq/URERKTtmUzw+9/DmDGQlASzZhkdkdwsu91e99v8M2fOcObMGSorK136jDbbJuzo0aMEBgY2qNIOHz4cgGPHjjUrwXWVvXv3MmrUqHZ7noiIiLhGSAjExMDrr8PMmdCtm9ERyc1orMB48OBB7r//fpc9o812USgqKqJ3794N2mvbioqK2urRIiIi4mFWr4YPPlByK83jlgc9iIiIiFztrrucL5HmaLMKbp8+fRqt0hYXF9f1i4iIiIi4WpsluCNGjOD48eNUV1fXaz9y5AgAw4YNa6tHi4iIiEgn1mYJbkREBOXl5aSnp9dr37x5M/7+/owZM6atHi0iIiIinVir1uDu3r2biooKysrKAOeOCLWJbHh4OD4+PkyePJlJkyYxe/ZsSktLCQgIwG6388knn5CcnNzk4RAiIiIiIjejVQnunDlzyM/PB5ynmKWlpZGWlobJZCIvL4/+/fsDkJGRwaJFi1iyZAnFxcUEBgaSkpLCk08+6bpPICIiIiJylVYluHl5ec26ztfXlw0bNrBhw4bWPEZERESkSe++6/zzt781Ng5xP222BldERESkLR09CvPnQ2Gh0ZGIu1GCKyIiIh3S0qXOo3yXLzc6EnE3SnBFRESkQ7Jancf3vvsuHD9udDTiTpTgioiISIf1u99B//4wb57RkYg7UYIrIiIiHVbXrmCzwc6dsHev0dGIu1CCKyIiIh3a9Okwbhy88gpcvmx0NOIOlOCKiIhIh2Yywe9/D//v/8EPPxgdjbiDVu2DKyIiIuJOQkLgxAnw8jI6EnEHSnBFRETEI3h5QWFhIQkJNrKzc6iq8sLb+zIhIUHYbAlYrVajQ5R2ogRXREREPEJBQQEPPPAUublvAjbABFSTk5PNgQNRZGamKsntJLQGV0RERDzC/Plra5LbsTiTW3CmOmPJzV1FQoLtuvcXFhYSFzePoUPDGTJkKkOHhhMXN49CHZXW4aiCKyIiIh4hOzsHZ+W2MWPYv38lZ87AHXc4v5h2NVV/PYsquCIiIuIRqqq8uFK5vZaZ/Hwv7rwTLBYYNQqio53H/R46dPPVX3EvquCKiIiIR/D2vgw4aDzJrSYg4DLr18PJk84dF06ehM8/h0GDblz9zc5e2WZxi+spwRURERGPEBISRE5OFs4q7LWyeOihIKZMadjjcMDKldev/jqrw9JRaImCiIiIeASbLYGAgIVAJlBd01oNZBIQsAibLaHR+0ymq6u/jamu6ZeOQhVcERER8QhWq5XMzNSafXBXXrMP7vW/JHaj6m9ISFCbxS2upwRXREREPIbVaiUpaW2L77PZEjhwIIrc3FXAGJy/5K4Gsmqqv6kujlTakhJcERER6fQaq/5evHiZ06eDeOMNbRHW0SjBFREREaFh9dfhgF//GlauhBkzwFtZU4ehL5mJiIiINMJkgnXrICcHkpKMjkZaQgmuiIiISBPuvx9iYmDxYigvNzoaaS63TnAzMzMxm82sWrWqru3YsWOMHz+enj17EhAQwPvvv29ghCIiIuLpVq2C8+dhbcu/uyYGcdsEt7q6mvj4eEJDQzFddWD0M888Q1hYGOfPnyc9PZ34+Hi+/fZbAyMVERERT3b33TB3rnO5wk8/GR2NNIfbLpd+7733GDduHMXFxTgcVzZePn78OKmpzq06goODCQwM5MSJE9x3331GhSoiIiIe7rXX4McfobLS6EikOdyygltUVMTGjRtZunRpg76wsDA++OADqqqqyMrK4vTp04SGhhoQpYiIiHQWPXrA1q3Oaq64P7dMcF977TVeeeUVunfvDlBviYLNZiMpKQkfHx8efPBBVqxYgZ+fn1GhioiIiIibMTzBTU5OxmKxYLFYCA8P5+uvv+abb77h2WefBcDhcNQtUaioqODRRx9l3bp1XLp0iZMnT7JhwwZ27Nhh5EcQERERETfS4gS3vLychIQEwsLCsFqtmM1mli9f3uS1c+fOxd/fHx8fH4KDg+vWz9aKiYmhrKyMsrIydu7cyRdffEFOTg5+fn5YrVZSU1NZvXo1sbGxHD9+nKqqKmbMmIHJZGLgwIFMmTKFPXv2tO7Ti4iIiIjHaXGCe+7cOTZt2kRlZSURERFA/SUEV4uMjGTLli0sW7aMjz/+mNGjRxMdHY3dbm/y/WfNmsWpU6c4fPgwhw4dYurUqbzwwgusX7+ee+65h9LSUjIyMnA4HOTn5/OXv/yFESNGtPRjiIiIiIiHavEuCgMGDKCkpARwfhksMTGx0et27drFvn37sNvtREVFAfDwww+Tn5/PvHnziIqKwmxumF/7+vri6+tb93O3bt3o3r07vXr1AiAlJYXXX3+d2NhYLBYLM2fOZNasWS39GCIiIiLioW5qm7Crt++61rZt27BYLMyYMaNee1xcHDNnziQrK6tZux8kXXM23uOPP87jjz/euoBFREREXKSkBL7+GiZONDoSuVab7YN79OhRAgMDG1Rphw8fDjhPJGvP7b0mTZpEly5d8Pf3x9/fH4Do6Giio6PbLQYRERHxHGvXwsaN8N13cNttRkfTcdjt9rrlqmfOnOHMmTNUuniD4TZLcIuKihg0aFCD9t69e9f1t6e9e/cyatSodn2miIiIeK558+Ddd2HZMuef0jyNFRgPHjzI/fff77JnGL5NmIiIiEhH1KsXLF4MmzZBTo7R0cjV2izB7dOnT6NV2uLi4rp+ERERkY5szhwYMADmzzc6ErlamyW4I0aM4Pjx41RXV9drP3LkCADDhg1rq0eLiIiItIuuXWHNGtixA/76V6OjkVptluBGRERQXl5Oenp6vfbNmzfj7+/PmDFj2urRIiIiIu3miSdg7Fh49VW4pq4nBmnVl8x2795NRUUFZWVlgHNHhNpENjw8HB8fHyZPnsykSZOYPXs2paWlBAQEYLfb+eSTT0hOTm7ycAgRERGRjsRkgnXr4MEHCwkLs/G//5tDVZUX3t6XCQkJwmZLwGq1Gh1mp9KqBHfOnDnk5+cDzlPM0tLSSEtLw2QykZeXR//+/QHIyMhg0aJFLFmyhOLiYgIDA0lJSeHJJ5903ScQERERMdi99xbg7/8Un376JmADTEA1OTnZHDgQRWZmqpLcdtSqBDcvL69Z1/n6+rJhwwY2bNjQmseIiIiIdAjz56/lzJk3gbFXtZqBseTmriIhwUZS0tom7y8sLCQhwUZ2tqq/rtBm++CKiIiIdBbZ2Tk4K7eNGUN29som7y0oKOCBB54iN1fVX1dRgisiIiJyk6qqvHAmpo0xU1DgxdixMHgwDBly5c9773VWf53Jbeuqv9KQElwRERGRm+TtfRlw0HiSW42Pz2UCA+HECdi5E2qOBai5t/XVX2mcTjITERERuUkhIUFAVhO9WTz6aBBJSfDll1BUBOfOwd//Dn/8I1gs16/+OqvD0hJKcEVERERuks2WQEDAQiATqN0MtxrIJCBgETZbQr3r+/SBBx6AuDi4/fba6m9jqmuqw9ISWqIgIiIicpOsViuZmak1OyGsvGYnhOt/SSwkJIicnCzqr8GtlVVTHZaWUIIrIiIi4gJWq7VVXwaz2RI4cCCK3NxVwBicv2CvBrJqqr+pLo7U8ynBFRERETHQzVR/pXFKcEVEREQM1lj19/Bh+NOf4KWXDAqqA9OXzERERETc0IED8PLLcOSI0ZF0PEpwRURERNzQc8/BPffAggVGR9LxKMEVERERcUO33AJvvgm7dsFf/2p0NB2LElwRERERN/XEExASAgkJUF194+vFSQmuiIiIiJsymWDtWvj6a/jwQ6Oj6TiU4IqIiIi4sfHjYepUWLgQfvnF6Gg6BiW4IiIiIm5u9Wr46Sf44gujI+kYtA+uiIiIiJsLCoLTp0FnPjSPKrgiIiIiHYCS2+ZTgisiIiIiHkUJroiIiIh4FCW4IiIiIuJRlOCKiIiIiEdRgisiIiLSQV2+bHQE7sktE9wJEybg4+ODxWLBYrEwceLEur6nnnqK2267je7duzNkyBA2bdpkYKQiIiIixnjnHXjwQR3h2xi3THBNJhPvv/8+ZWVllJWVsW/fvrq+pUuX8uOPP1JaWsrWrVt58cUXycvLMzBaERERkfY3YgT8939DSorRkbgft0xwARwOR6PtgYGBeHs7z6fw8vKie/fuWCyW9gxNRERExHAPPQT/+q86wrcxbpvgxsfH4+fnx8SJEzl8+HC9vpiYGHx8fBg3bhx/+MMf6Nu3r0FRioiIiBhnzRr48UfncgW5wi0TXJvNxvfff8/p06cJDw/nN7/5DaWlpXX9ycnJVFRUYLfbiYuL44cffjAwWhERERFj3HcfzJoFK1dCSYnR0bgPwxPc5OTkui+ThYeHAzB69Gi6detG165diY+Pp2/fvnz55Zf17jObzUybNo1x48axfft2I0IXERERMdyyZfDLL4WMHz+PoUPDGTJkKkOHhhMXN4/CwkKjwzNEixPc8vJyEhISCAsLw2q1YjabWb58eZPXzp07F39/f3x8fAgODiY1NbXeNTExMXVfJtu5c2eLP0BVVRW33npri+8TERER8QRmcwFdu0Zx9Oh0cnJ2cPLkdnJyPmLz5umEhkZ1yiS3xQnuuXPn2LRpE5WVlURERADOXQ8aExkZyZYtW1i2bBkff/wxo0ePJjo6Grvd3uT7X7hwgb179/LLL79w6dIlNm7cSEFBAaGhoZw9e5b09HQqKiqoqqriww8/JCsri0mTJrX0Y4iIiIh4hPnz11JS8iYwFqjNyczAWHJzV5GQYLvu/YWFhcTFeVb117ulNwwYMICSmkUeRUVFJCYmNnrdrl272LdvH3a7naioKAAefvhh8vPzmTdvHlFRUZjNDfPryspKFi5cyIkTJ+jSpQsjR45k165d9OjRg59//pmNGzcya9YszGYzw4YNY/v27fj7+7f0Y4iIiIh4hOzsHKCpJHYM2dkrm7y3oKCABx54itzcN2vewwRUk5OTzYEDUWRmpmK1Wl0fdBtrcYJ7taa28gLYtm0bFouFGTNm1GuPi4tj5syZZGVlERoa2uC+vn378o9//KPR9+zXrx9/+9vfWhXrpEmT6NKlC/7+/nUJcXR0NNHR0a16PxERERF3UFXlxZXK7bXMNf2wahXs2wdDhjhfgwfDH/+4tia5HVvvnqurv0lJa10ar91ur/tt/pkzZzhz5gyVlZUufcZNJbjXc/ToUQIDAxtUaYcPHw7AsWPHGk1w28revXsZNWpUuz1PREREpD14e18GHDSe5FbX9MOAAeDn5zwcYssWuHgRoPXV39ZqrMB48OBB7r//fpc9o812USgqKqJ3794N2mvbioqK2urRIiIiIp1GSEgQkNVEb1ZNP8TEQGoqHDoE5eWQnw/+/s2r/nY0hm8TJiIiIiKtZ7MlEBCwEMgEqmtaq4FMAgIWYbMlNLjHbIb+/aFHj9rqb2OuVH87mjZbotCnT59Gq7TFxcV1/SIiIiJyc6xWK5mZqSQk2MjOXklVlRfe3pcJCQnCZrv+l8RCQoLIycmi/hrcWleqvx1NmyW4I0aMwG63U11dXW8d7pEjRwAYNmxYWz1aREREpFOxWq2t+jKYzZbAgQNR5OauAsbg/OV+NZBVU/1Nvf4buKk2W6IQERFBeXk56enp9do3b96Mv78/Y8aMaatHi4iIiEgz1FZ/Y2MzCAqawuDBUwkKmkJsbEaH3SIMWlnB3b17NxUVFZSVlQHOHRFqE9nw8HB8fHyYPHkykyZNYvbs2ZSWlhIQEIDdbueTTz4hOTm5ycMhRERERKT9tLb6685aleDOmTOH/Px8wHmKWVpaGmlpaZhMJvLy8ujfvz8AGRkZLFq0iCVLllBcXExgYCApKSk8+eSTrvsEIiIiIuJyDgdcuAA9exodScu1KsHNy8tr1nW+vr5s2LCBDRs2tOYxIiIiImKQl1+GL76A7GzoaL941zZhIiIiItLAtGnw1VdwzdepOgQluCIiIiLSwMMPw7/8CyxcCC4+SbfNKcEVERERkUatXg25uZCYaHQkLaMEV0REREQaNWIEPP00LF/uPN63o1CCKyIiIiJNWrECSkpg/XqjI2k+JbgiIiIi0qS774YXXgCbDQoLjY6medrsqF4RERER8QwLFzr3xTV3kNKoElwRERERua4+feD3vzc6iubrIHm4iIiIiEjzKMEVEREREY+iBFdEREREPIoSXBERERHxKEpwRURERMSjKMEVERERkRb77jv48kujo2ictgkTERERkRZ75RXIyXG+unQxOpr6VMEVERERkRZ74w3IzYXERKMjaUgJroiIiIi02IgR8PTTsHw5lJcbHU19WqIgIiIiIq2yYgXY7YU89JCNS5dyqKrywtv7MiEhQdhsCVitVkPiUoIrIiIiIq3i41NAt25PcejQm4ANMAHV5ORkc+BAFJmZqYYkuVqiICIiIiKtMn/+Wi5ceBMYizO5BWd6OZbc3FUkJNiue39hYSFxcfN44okXXRqXElwRERERaZXs7BxgTBO9Y2r6G1dQUEBoaBSbN08nL+9tl8alJQoiIiIi0ipVVV5cqdxey1zTf8U778Add8CQIbBmzVpyc2urvwddGpdbVnAnTJiAj48PFosFi8XCxIkT6/qOHTvG+PHj6dmzJwEBAbz//vsGRioiIiLSeXl7XwYcTfRW1/Q7VVbC4sUQGQlDh8IHH1yv+ntz3DLBNZlMvP/++5SVlVFWVsa+ffvq+p555hnCwsI4f/486enpxMfH8+233xoYrYiIiEjnFBISBGQ10ZtV0+/UpQsUFcH//R98/jn063e96u/NccsEF8DhaPxfA8ePHyc6OhqA4OBgAgMDOXHiRHuGJiIiIiKAzZZAQMBCIBOormmtBjIJCFiEzZZQ73qTCfr1g/HjoU+f61V/b47bJrjx8fH4+fkxceJEDh8+XNceFhbGBx98QFVVFVlZWZw+fZrQ0FADI/Vcdrvd6BA6HI1Z62jcWk5j1joat5bTmLVOZxk3q9VKZmYqsbEZBAVNYfDgqQQFTSE2NuOGW4Rdv/p7c9wywbXZbHz//fecPn2a8PBwfvOb31BaWlrXl5SUhI+PDw8++CArVqzAz8/P4Ig9U2f5y+lKGrPW0bi1nMasdTRuLacxa53ONG5Wq5WkpLUcO7aTEye2c+zYTpKS1t5w/9vGq7+uYXiCm5ycXPdlsvDwcABGjx5Nt27d6Nq1K/Hx8fTt25fMzEwqKip49NFHWbduHZcuXeLkyZNs2LCBHTt2GPwpRERERKQlrq7+Dhw416Xv3eIEt7y8nISEBMLCwrBarZjNZpYvX97ktXPnzsXf3x8fHx+Cg4NJTU2td01MTEzdl8l27tzZ5HMdDgc5OTlUVVUxY8YMTCYTAwcOZMqUKezZs6elH0NEREREDFZb/U1P3+jS921xgnvu3Dk2bdpEZWUlERERgHPXg8ZERkayZcsWli1bxscff8zo0aOJjo6+btn+woUL7N27l19++YVLly6xcePGmo2AQwkICKC0tJSMjAwcDgf5+fn85S9/YcSIES39GCIiIiLioVp80MOAAQMoKSkBoKioiMTExEav27VrF/v27cNutxMVFQXAww8/TH5+PvPmzSMqKgqzuWF+XVlZycKFCzlx4gRdunRh5MiR7Nq1ix49egCQkpLC66+/TmxsLBaLhZkzZzJr1qwm47148SLg3H1BWubChQscPOjajZc9ncasdTRuLacxax2NW8tpzFpH49YytXlabd520xw3obCw0GEymRzLly9v0Ddr1ixH9+7dHZcvX67XbrfbHSaTyfHll1/ezKObbevWrQ6ce1DopZdeeumll1566eXGr61bt7ok/2uzo3qPHj1KYGBggyrt8OHDAeeJZO2xvddjjz3G1q1bGTBgAD4+Pm3+PBERERFpmYsXL/L999/z2GOPueT92izBLSoqYtCgQQ3ae/fuXdffHvr27UtMTEy7PEtEREREWmfcuHEuey/DtwkTEREREXGlNktw+/Tp02iVtri4uK5fRERERMTV2izBHTFiBMePH6e6uv7JFEeOHAFg2LBhbfVoEREREenE2izBjYiIoLy8nPT09Hrtmzdvxt/fnzFjxrTVo0VERESkE2tVgrt7927S09P56KOPAOeOCOnp6aSnp9ftXzZ58mQmTZrE7NmzSUxMZP/+/Tz33HN88skn2Gy2Jg+HcJXmnKIm9X322WeYzeZGX9nZ2UaHZ7iWnOJ38OBBJk6ciMVioVevXkyfPp28vLx2jtg9NHfcYmNjG517QUFBBkRtrE8//ZRnnnmGwYMH4+vry5133sm0adMa3VNTc82puWOmeVbfoUOHCA8P5+6776Zbt2706dOHBx54gOTk5AbXaq45NXfMNNeuLzExEbPZjMViadDnirnWql0U5syZQ35+PuA8xSwtLY20tDRMJhN5eXn0798fgIyMDBYtWsSSJUsoLi4mMDCQlJQUnnzyydY8tkUiIyP56quv+I//+A8GDx5McnIy0dHRVFdXEx0d3ebP78hWr17NI488Uq9t6NChBkXjPmpP8Rs5ciQREREkJiY2+g+1b7/9lgkTJjBq1CjS0tK4ePEiS5Ys4aGHHuLQoUP07dvXgOiN09xxA/Dx8WH//v0N2jqb9957j8LCQuLj4xk6dCiFhYW89dZbjB07lj179tT9/dRcu6K5YwaaZ1e7cOEC/fv3JyYmBn9/f8rLy0lOTubpp5/m+++/Z9GiRYDm2tWaO2agudaUM2fO8Oqrr3LHHXdQWlpar89lc80lu+m6mZ07dzpMJpMjJSWlXntYWJjD39+/weET4rR//36HyWRy/PnPfzY6FLd37ty5Jg85mTFjhsPPz89RVlZW15afn++45ZZbHPPnz2/PMN3O9cbtmWeecVgsFgOicj9nz55t0FZeXu647bbbHBMnTqxr01y7orljpnnWPGPHjnX079+/7mfNtRu7dsw015r2+OOPO6ZNm+aIjY113HrrrfX6XDXXPHKbsG3btmGxWJgxY0a99ri4OH766SeysrIMiqxjcDgcRofg9poao6qqKnbs2MH06dO59dZb69r79+/PI488wrZt29orRLd0o7mluefk5+fXoM3X15fAwEB+/PFHQHPtWs0Zs1qaZzfWp08fvL2dv+TVXGueq8esluZaQ1u3buXAgQO88847DcbHlXPNIxPc5pyiJk17/vnn6dKlCz169GDy5Mn8/e9/NzqkDiM3N5eff/6ZESNGNOgbPnw43333HZcuXTIgso7h4sWL3H777Xh7e3PXXXfxu9/9jpKSEqPDcgu159rXLhfSXLuxa8esluZZQw6Hg6qqKgoLC/mv//ov9uzZw6uvvgporjXlemNWS3OtvrNnzzJ37lzWrFnDHXfc0aDflXOtzU4yM5K7nKLW0fTs2ZO5c+cyYcIE+vTpw6lTp1i7di0TJkxg586dhIWFGR2i26udW7Vz7Wq9e/fG4XBQUlJCv3792js0tzdy5EiCg4PrthD87LPPWL9+PZ9++in/+Mc/8PX1NThCYz3//PNcvHixbn2f5tqNXTtmoHnWlNmzZ/OHP/wBAC8vL9atW8fs2bMBzbWmXG/MQHOtMc8//zxBQUH89re/bbTflXPNIxNcaZ2RI0cycuTIup/HjRtHREQEw4cPZ/78+UpwpU3NnTu33s+PPvoowcHBPPHEEyQmJvLSSy8ZFJnxFi9ezJ/+9Cf+8z//k+DgYKPD6RCaGjPNs8YtWrSI5557joKCArZv387LL7/Mzz//zPz5840OzW3daMw01+pLT09nx44dHD58uF2e55EJrk5Rc50ePXoQHh7Oe++9xy+//ELXrl2NDsmt1c6t2rl2teLiYkwmE7169WrvsDqsiIgIfH19O/W6+eXLl7Nq1SrefPNN5syZU9euuda0psasKZpncNddd3HXXXcBzm0+wfmPhGefYdKn9QAAA5xJREFUfVZzrQlNjdm///u/Y7VaG72ns8618vJyXnjhBV588UX69evH+fPnAeqWG1y4cAFvb2+XzjWPXIOrU9TaRlvvXewJAgIC8PHx4X/+538a9B05coR7772XW265xYDIOiaHw9Hg73Fnsnz58rrXggUL6vVprjXuemPWlM4+zxozevRoqqqq+Oc//6m51ky1Y3a9/Vo761w7d+4cBQUFrFu3jt69e9e9UlJSqKiooFevXjz99NMMGjTIZXPNIxNcnaLmOiUlJXz00UcEBwfrf2DN4O3tzZQpU8jIyKC8vLyu/YcffmD//v1ERkYaGF3HU3t4TGhoqNGhtLs33niD5cuXs3jxYhYvXtygX3OtoRuNWVM68zxryv79+/Hy8iIgIEBzrZmuHrOmdNa5dvvtt7N//34+++yzutf+/ft57LHH+NWvfsVnn33GypUr8fLyctlcMzk8dA+Lxx57rO6gh4CAAOx2O4mJiXUHPkhDMTExDBw4kFGjRtG7d29OnTrFW2+9RV5eHrt37+bXv/610SEabvfu3VRUVFBWVsazzz7LjBkz6rajCw8Px8fHhxMnTjB69GhGjRrFggUL6japPn/+PIcOHeqUS2RuNG4FBQX827/9GzNnzuSee+7B4XDw+eef8/bbbzNo0CCysrI61ebob731FvPmzWPy5MksXbq0wVY6Y8eOBdBcu0pzxiw/P1/z7BrPPfccPXr0YPTo0fTr149z586RlpbGhx9+SEJCAmvWrAE0167WnDHTXGue2NhY/vznP1NWVlbX5rK51vLteTuG8vJyx0svveS4/fbbHV27dnWMHDnSkZqaanRYbm3NmjWO4OBgR8+ePR3e3t4OPz8/x/Tp0x1fffWV0aG5jQEDBjhMJpPDZDI5zGZzvf/Oz8+vu+7rr792TJw40eHr6+vo0aOHIzIy0vHPf/7TwMiNdaNxKykpcURGRjoGDhzo6Natm6Nr166OIUOGOBYsWOAoLS01Ovx2N2HChHrjdPXLbDbXu1Zzzak5Y6Z51lBSUpJj/PjxDqvV6ujSpYujV69ejkceecSRnJzc4FrNNafmjJnmWvPExsY2ehiGK+aax1ZwRURERKRz8sg1uCIiIiLSeSnBFRERERGPogRXRERERDyKElwRERER8ShKcEVERETEoyjBFRERERGPogRXRERERDyKElwRERER8ShKcEVERETEoyjBFRERERGPogRXRERERDzK/wfzS9AZDasJOgAAAABJRU5ErkJggg==",
      "text/plain": [
       "PyPlot.Figure(PyObject <matplotlib.figure.Figure object at 0x7f784abd15d0>)"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "1-element Array{Any,1}:\n",
       " PyObject <matplotlib.lines.Line2D object at 0x7f784aaf2e50>"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "semilogy(Q, derr, \"o--\")"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Julia 0.4.5",
   "language": "julia",
   "name": "julia-0.4"
  },
  "language_info": {
   "file_extension": ".jl",
   "mimetype": "application/julia",
   "name": "julia",
   "version": "0.4.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}
